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ABSTRACT 

The fraction of stellar mass contained in globular clusters (GCs), also measured by number as the 
specific frequency, is a fundamental quantity that reflects both a galaxy's early star formation and 
its entire merging history. We present specific frequencies, luminosities, and mass fractions for the 
globular cluster systems of 100 early-type galaxies in the ACS Virgo Cluster Survey. This catalog 
represents the largest homogeneous catalog of GC number and mass fractions across a wide range of 
galaxy luminosity (—22 < Mb < —15). We find that 1) GC mass fractions can be high in both giants 
and dwarfs, but are universally low in galaxies with intermediate luminosities (—20 < Mb < —17). 
2) The fraction of red GCs increases with galaxy luminosity, but stays constant or decreases for 
galaxies brighter than M z = —22. As a result, although specific frequencies for blue and red GCs are 
both higher in massive galaxies, the behavior of specific frequency across galaxy mass is dominated 
by the blue GCs. 3) The GC fractions of low-mass galaxies exhibit a dependence on environment, 
where dwarf galaxies closer to the cluster center have higher GC fractions. Nearly all dwarfs with 
high GC fractions are within 1 Mpc of the cD galaxy M87, presenting the first strong evidence that 
GC formation in dwarf galaxies is biased toward dense environments. 4) GC formation in central 
dwarfs is biased because their stars form earliest and most intensely. Comparisons to early-type 
dwarf galaxies in the Millennium Simulation show that central dwarfs are likely to have older stellar 
populations and form more of their stars at higher star formation rates (SFRs) and star formation 
rate surface densities. In addition, the SFR surface density in simulated dwarfs peaks before the total 
SFR, naturally producing GC populations that are older and more metal-poor than the field stars. 5) 
Dwarfs within ~ 40 kpc of the giant ellipticals M87 and M49 are red for their luminosities and have 
few or no GCs, suggesting that they have been tidally stripped and have contributed their GCs to the 
halos of their giant neighbors. The central dwarfs with high GC mass fractions are thus likely to be 
the survivors most similar to the protogalaxies that assembled the rich M87 globular cluster system. 
Subject headings: galaxies: elliptical and lenticular, cD — galaxies: dwarf — galaxies: halos - 
galaxies: evolution — galaxies: star clusters - globular clusters: general 
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1. INTRODUCTION 

Globular clusters (GCs) constitute a small fraction of 
the stellar mass in galaxies, but their ubiquity, relative 
simplicity, and old ages make them the most prominent 
representatives of a bygone epoch of galaxy formation. 
GCs are made of stars that are among the oldest in 
galaxies, and they can be observed at large distances 
(e.g., Blakeslee et al. 2003a). These old star clusters are 
thus unique, both intrinsically and observationally, for 
understanding the early, intense star-forming episodes 
that mark galaxy formation. 

In the local universe, we see massive star clusters form- 
ing wherever there are high star formation rate sur- 
face densities (Larsen & Richtler 2000), providing a 
connection that suggests the properties of star cluster 
populations — age, metallicity, mass — should scale quite 
closely with field stars formed in the same events. The 
properties of globular cluster systems do in fact correlate 
strongly with the properties of the field stars of their host 
galaxies. The mean metallicities of GC systems have 
long been known to scale with the metallicity of their 
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host (van den Bergh 1975; Brodie & Huchra 1991), and 
the mean metallicities of both the metal-rich and metal- 
poor subpopulations also correlate with the luminosity 
and mass of the host galaxy (Larsen et al. 2001; Peng 
et al. 2006a and references therein). However, if GC sys- 
tems directly followed the underlying field light in every 
way, they might be less interesting. For instance, al- 
though the metallicities of GC systems may track those 
of galaxies, they are consistently offset to lower values by 
~ 0.5-0.8 dex in [Fe/H] (e.g., Jordan et al. 2004a, Lotz 
et al. 2004). Most conspicuously, even the most massive 
and metal-rich galaxies have GC systems dominated by 
metal-poor star clusters ([Fe/H]< —1). This suggests a 
disconnect between the formation of "halo" stellar pop- 
ulations and the bulk of the galaxy. 

One of the most studied aspects of this GC-galaxy du- 
ality concerns the specific frequency of globular clusters, 
or the number of GCs per unit stellar luminosity. Spe- 
cific frequency, Sn, was introduced by Harris & van den 
Bergh (1981) and is defined as the number of GCs nor- 
malized to a galaxy luminosity of My = — 15. The pur- 
pose of studying Sn across galaxies of different masses, 
morphologies, and environments is, in the words of that 
initial paper, "to investigate whether there is in fact a 
'universal' and uniform capability for globular cluster for- 
mation" . This simple quantity, and similar ones related 
to it, turn out to be extremely interesting galaxy diag- 
nostics. It appears true that for galaxies above a certain 
mass there is a nearly universal capability to form glob- 
ular clusters, but this process is not uniform across all 
galaxies, at least as seen in comparison to the field stars. 
Specific frequencies of spiral galaxies like the Milky Way 
are generally 0.5-1 (Goudfrooij et al. 2003; Rhode & Zepf 
2004; Chandar, Whitmore & Lee 2004), although they 
have a mean of ~ 4 when normalized only to bulge lu- 
minosity (Cote et al. 2000). The specific frequencies of 
massive ellipticals are 2-6, and those of some cD galaxies 
such as M87 (N4486) can be well in excess of 10. These 
trends are apparent even when the number of GCs is nor- 
malized to stellar mass as opposed to stellar luminosity 
(Rhode, Zepf & Santos 2005). Dwarf elliptical galaxies 
(dEs), whose GC systems are predominantly metal-poor, 
can also have high Sn similar to those of giant ellipticals 
(Durrell et al. 1996; Miller et al. 1998; Lotz et al. 2004; 
Miller & Lotz 2007; Puzia & Sharina 2007), as can some 
dwarf irregulars (Seth et al. 2004), suggesting the possi- 
bility that the halos of large galaxies were formed mainly 
through the accretion of dwarf-like objects (Searle & Zinn 
1978; Cote, Marzke & West 1998; Cote et al. 2000). 

A central question in the study of GC systems is: How 
do we understand different GC fractions in the context of 
galaxy assembly? The formation of globular cluster sys- 
tems has been particularly tied to the formation of mas- 
sive elliptical galaxies, in which GCs are often present in 
large numbers and where GCs are most easily observed. 
The mergers and accretion events expected during the hi- 
erarchical assembly of these galaxies must also be able to 
form their GC systems. Observations and simulations of 
elliptical galaxy formation are creating a picture in which 
the stars form early and quickly (e.g., Kodama et al. 
1998), mimicking the traditional "monolithic collapse" 
scenario, but where the assembly of these stars into a 
single galaxy continues until late times through largely 
dissipationless mergers (De Lucia et al. 2006; De Lucia 



& Blaizot 2007), and star formation at late times is sup- 
pressed by energy feedback (Springel et al. 2005; Croton 
et al. 2006). 

Intertwined with the issue of galaxy formation is that 
of the formation efficiency of the GCs themselves: Why 
do globular clusters form with different efficiencies with 
respect to their light in different galaxies? Blakeslee, 
Tonry & Metzger (1997) and Blakeslee (1999) studied 
the GC systems of brightest cluster galaxies (BCGs) in 
Abell galaxy clusters and found that the number of GCs 
scaled with the velocity dispersion of the galaxy clus- 
ter rather than with the luminosity of the BCGs, sug- 
gesting that GC formation is closely linked to the total 
mass of the system, i.e., Sn M./L. In a similar vein, 
McLaughlin (1999a) examined the high Sn in M87 and 
found that the large number of GCs it possessed was not 
anomalous when normalized to the total baryonic mass 
(including the hot X-ray gas) rather than just to the stel- 
lar mass. McLaughlin (1999a) defined a "universal" GC 
formation efficiency of e = 0.26%, where e is the fraction 
of the baryonic mass that ends up in globular clusters. 
Kravtsov & Gnedin (2005) studied the formation of the 
GC system in a high resolution hydrodynamic simula- 
tion and also found that the mass in GCs was directly 
proportional to the total halo mass of the galaxy. If to- 
tal mass drives GC formation, then it is the variation in 
converting baryons into field stars that drives trends in 
specific frequency. 

The connection between GCs and galactic mass (bary- 
onic or total) makes it tempting to try to explain them 
using simulations of dark matter and galaxy assembly. 
Beasley et al. (2002) simulated the color distributions 
of GCs using semi-analytic models of galaxy forma- 
tion, while the aforementioned simulations of Kravtsov 
& Gnedin (2005) formed a Milky Way mass galaxy in 
detail. Moore et al. (2006) identified metal-poor GCs 
with early collapsing dark matter peaks in cosmological 
N-body simulations. An empirical connection between 
the total mass of a system and the mass contained in its 
globular clusters does not, however, explain why the star 
formation histories of the GCs and the field should be 
different. 

One clue is that in the local Universe, high mass frac- 
tions in massive star clusters occur in regions of high 
star formation surface density (Larsen & Richtlcr 2000). 
Thus, one explanation for high-Sjy galaxies might be that 
they formed more of their stars in high efficiency events. 
Given that the star formation rate (SFR) in galaxies was 
most intense at very early times (z > 2), variations in 
Sn could result from different times of formation, espe- 
cially when coupled with a sharp star formation cutoff 
at reionization (Santos 2003). In this scenario, with Sn 
a function of formation time, GC formation is biased 
towards the earliest collapsing halos which can create a 
large fraction of their stars at high efficiency before reion- 
ization. Low mass halos in dense environments collapse 
earlier (Gao, Springel & White 2005; Diemand, Madau & 
Moore 2005), and these are also the ones most susceptible 
to heating through photoionization and feedback. One 
expectation of this scenario is that Sn will be "biased" 
towards dense environments (West 1993), especially in 
low mass galaxies. The simulations of Moore et al. (2006) 
suggested that GCs and satellites would be more highly 
clustered in dense regions for higher rcdshifts of reion- 
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ization, and the same idea can be applied to Sn at fixed 
reionization rcdshift. 

Previous observational studies have shown some ev- 
idence of an environmental dependence for Sn- West 
(1993) showed that the mean Sn of elliptical galaxies 
correlates with the local galaxy density. Blakeslee (1997, 
1999) found that Sn in brightest cluster galaxies (BCGs) 
scaled with properties that reflect cluster density (cluster 
velocity dispersion, X-ray temperature, and X-ray lumi- 
nosity), and that cluster galaxies closer to their cluster's 
X-ray center have higher Sn- West et al. (1995) also 
found a correlation between the Sn of BCGs and the X- 
ray luminosity of the cluster, interpreting it as evidence 
for a population of intcrgalactic globular clusters. On the 
other hand, Lotz et al. (2004) and Miller & Lotz (2007) 
did not find an obvious correlation between Sn and clus- 
tercentric radius for dEs in the Virgo and Fornax clus- 
ters, although the number of GCs in their galaxies were 
small and Sn errors were quite large. Also, Spitler et al. 
(2007) recently found that the relatively isolated ellipti- 
cal NGC 821 has an Sn comparable to cluster galaxies 
of like luminosities. 

The specific frequencies (or formation efficiencies) of 
globular cluster systems is clearly a fundamental prop- 
erty of galaxies. However, accurately measuring Sn is 
traditionally fraught with uncertainty despite the fact 
that at its most basic level it requires only simple count- 
ing of GCs. In practice, observations need to be deep 
enough to observe past the mean of the GC luminosity 
function, and GC selection needs to be efficient enough 
so that contaminants do not overwhelm the GCs in lower 
mass galaxies. The total magnitude and distance of the 
galaxy also needs to be known to establish its luminosity 
(or mass). Studies of larger galaxies also benefit from 
wide-field coverage so that enough of the GC system 
is sampled to minimize extrapolation errors. It is be- 
cause of these limitations that homogeneous surveys of 
Sn across a wide range of galaxy mass are difficult to 
conduct. To put together a complete picture of globu- 
lar cluster system and galaxy formation, it is important 
to study galaxies at all masses in the same way — at low 
masses we study the survivors of hierarchical assembly, 
and at high masses we study its final products. 

With this in mind, we have undertaken a careful study 
of the formation efficiencies for the GC systems of 100 
early-type galaxies in the ACS Virgo Cluster Survey 
(ACSVCS; Cote et al. 2004). The deep, high-resolution, 
relatively wide-field imaging provides the most complete 
census to date of GCs in early-type galaxies over a large 
range in galaxy luminosity (—22 < Mb < —15). In this 
paper, we present the specific frequencies and other re- 
lated quantities of the ACSVCS galaxies and quantify 
their trends as a function of host galaxy properties. 

2. OBSERVATIONS AND DATA 

The ACS Virgo Cluster Survey (ACSVCS) is a large 
program to image 100 early-type galaxies in the Virgo 
Cluster with the HST Advanced Camera for Surveys 
(ACS; Ford et al. 1998). The survey is described in 
detail in Paper I (Cote et al. 2004), and the data re- 
duction techniques are outlined in Jordan et al. 2004b 
(Paper II), but we briefly summarize them here. We 
obtained images in the F475W [g) and F850LP (z) fil- 
ters of galaxies selected to be early-type and with con- 



firmed cluster membership in the Virgo Cluster Catalog 
(VCC) of Binggeli, Sandage, & Tammann (1985). Our 
sample excludes the Southern Extension, has Bt < 16, 
and the galaxies are morphologically classified as E, SO, 
dE, dE,N, dSO, or dS0,N. This early-type galaxy sam- 
ple is magnitude limited for the brightest 26 galaxies 
(Bt < 12.15 or Mb < —19.10), and contains a represen- 
tative sample of galaxies for Mb < —15.11. The color 
distributions of the GC systems were presented in Peng 
et al. (2006a; Paper IX), and their size distributions in 
Jordan et al. (2005; Paper X). We have measured dis- 
tances for 84 galaxies from the method of surface bright- 
ness fluctuations (Mei et al. 2005a, 2005b, 2007; Papers 
IV, V, and XIII). For galaxies where an SBF distance 
could not be measured, we adopt a distance to the Virgo 
Cluster of D = 16.5 Mpc with a distance modulus of 
31.09 ± 0.03 mag from Tonry et al. (2001), corrected by 
the final results of the Key Project distances (Freedman 
et al. 2001; see also discussion in Mei et al. 2005b). The 
surface brightness profiles, total magnitudes, and colors 
of the sample galaxies in g and z were described in Fer- 
rarese et al. (2006a; Paper VI), and the properties of 
the GC luminosity functions for 89 of our sample galax- 
ies were presented in Jordan et al. (2007a, 2007b, Pa- 
per XII). The galaxy nuclear properties were presented 
in Cote et al. (2006; Paper VIII) and their connection 
to supermassive black holes in Ferrarese et al. (2006b). 
These data have also been analyzed for ultra-compact 
dwarf galaxies (Hasegan et al. 2005; Paper VII), diffuse 
star clusters (Peng et al. 2006b; Paper XI), the connec- 
tion between GCs and low mass X-ray binaries (Jordan 
et al. 2004c, Sivakoff et al. 2007), and color-magnitude 
relations in GC systems (Micskc et al. 2006; Paper XIV). 
Together, these data create the best opportunity to date 
to study the formation efficiencies of GC systems in early- 
type galaxies, and this paper will refer often to the quan- 
tities measured in the preceding papers. 

2.1. Data Reduction and Control Fields 

Each galaxy was imaged with the Wide Field Chan- 
nel (WFC) of the ACS. We reduced the ACS/WFC im- 
ages using a dedicated pipeline described in Jordan et al. 
(2004b, Paper II; sec also Blakeslee et al. 2003b). We 
produced the science images by combining and cleaning 
them of cosmic rays using the Pyraf routine multidrizzle 
(Koekemoer et al. 2002). We then subtracted a model 
of the galaxy light and used the source detection pro- 
gram SExtractor (Bertin & Arnouts 1996) to detect and 
mask sources and remove residual background. Our fi- 
nal object detection includes estimates of both the image 
noise and noise due to surface brightness fluctuations — 
ignoring the latter results in many false detections in 
the bright central regions of the galaxy — and objects 
are only included in the final catalog if they are de- 
tected in both filters. After rejecting very bright or elon- 
gated objects to eliminate obvious foreground stars and 
background galaxies, and passing our catalog through 
a generous color cut, we use the program KINGPHOT 
(Jordan et al. 2005) to measure magnitudes and King 
model parameters. KINGPHOT fits King (1966) model 
surface brightness profiles convolved with the filter- and 
spatially-dependent point spread function (PSF). Mag- 
nitudes and colors are corrected for foreground extinc- 
tion using the reddening maps of Schlegel, Finkbeiner, & 
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Davis (1998) and extinction ratios for a G2 star (Paper 
II; Sirianni et al. 2005). For the purposes of this paper, 
whenever we refer to g and z, we mean the HST/ACS 
magnitudes 5475 and z 850 . Magnitudes in the Sloan Dig- 
ital Sky Survey (SDSS) system will be designated ex- 
plicitly, e.g., as in g s dss- The low redshifts of the galax- 
ies being studied (—575 to 2284 km s _1 ) means that K- 
corrections are negligible and thus we do not ^-correct 
any of the magnitudes presented in this paper. 

In any study of extragalactic star clusters, it is im- 
portant to quantify and correct for background contam- 
ination. To that end, we have also reduced 16 blank, 
high-latitude control fields taken from the ACS Pure Par- 
allel Program (GO-9488, PI: Ratnatunga; and GO-9575, 
PI: Sparks). For each galaxy, we have "customized" the 
control sample to mimic the spatially varying detection 
efficiency that is a function of the surface brightness of 
the unresolved galaxy light. Details of this procedure are 
given in Peng et al. (2006a, b; Papers IX and XI). 

We select probable globular clusters using their sizes 
and magnitudes. We can assign a probability that any 
object is a GC based on position in the r^-z diagram and 
the locus of contaminants in the same plane. Objects 
with probabilities greater than 0.5 are included in our 
GC sample, although the exact choice of the cutoff value 
does not significantly affect our results. We describe the 
details of this selection in Peng et al. (2006a, Paper IX) 
and Jordan et al. (2008, in preparation). 



2.2. WFPC2 Parallel Fields 

In addition to ACS/WFC imaging of our target galax- 
ies, we have also acquired parallel imaging of 100 
"blank" fields using the Wide Field Planetary Camera 
2 (WFPC2) in the F606W (wide V) and F814W (J) fil- 
ters. The separation between ACS and WFPC2 in the 
HST focal plane is approximately 5.'8, which is 29 kpc 
at the mean distance to the Virgo cluster (16.5 Mpc). 
For most of the ACSVCS galaxies, this distance puts 
the WFPC2 field well outside of the target galaxy, and 
thus one motivation for taking them is to search for in- 
tergalactic GCs (e.g., West et al. 1995, Jordan et al. 
2003, Williams et al. 2007). In this paper, we use the 
WFPC2 parallel images instead to constrain the total 
number of GCs in the larger, more luminous galaxies, 
where the GC system is still detectable 29 kpc from the 
center. The WFPC2 images were reduced using a mod- 
ified form of the PyRAF pipeline written by Alasdair 
Allan 18 . We implemented the cosmic ray cleaning algo- 
rithm LACOSMIC (van Dokkum 2001) and a geometric 
distortion correction (Bagget et al. 2002) in the PyRAF 
pipeline. Sources were selected from a SExtractor cata- 
log having magnitudes and colors of GCs. The sizes and 
magnitudes of GC candidates were also measured with 
KINGPHOT. To account for sample contamination, we 
selected 10 blank WFPC2 fields from the HST archive of 
similar depth and observed with the same filters as our 
WFPC2 data. We analysed them in the same manner 
as our data, using WFPC2 point spread functions pro- 
vided by P. B. Stetson. We find that there are on average 
5.8 ± 2.5 contaminant sources per WFPC2 field. 

18 http://www.astro.ex.ac.uk/pcoplc/aa/pagcs/computing/pyraf. 



3. CALCULATING SPECIFIC FREQUENCIES, 
LUMINOSITIES, AND MASSES 

3.1. Total Numbers of GCs 

The determination of specific frequency requires 
knowledge of the total numbers of GCs, the total ap- 
parent magnitude of the galaxy (traditionally in the V 
band), and the distance. For the distances, we use the 
SBF-dctcrmined distances presented in Mei et al. (2007; 
Paper XIII), using their polynomial calibration. The 
measurement of the total number of GCs is in principle 
a simple task of counting, but is not quite so straight- 
forward in practice. Our images are deep enough that 
we generally detect the brightest ~ 90% of the GC lumi- 
nosity function (GCLF), limiting the uncertainties from 
extrapolations in luminosity. Nevertheless, a full count 
of the GCs within the ACS field of view involves know- 
ing our level of completeness, which is a function of the 
GC magnitude, size, and the local background flux. We 
also need to know the form of the GCLF so that we can 
estimate the number of GCs missed. 

Fortunately, both the level of completeness and the 
GCLF can be determined for our data. We empirically 
derived the completeness of our data as a function of 
magnitude, background flux, and GC size, using simula- 
tions of nearly 5 million GCs placed in actual ACSVCS 
images for galaxies of different surface brightnesses. Us- 
ing these data, we can nonparametrically estimate the 
expected detection probability of any GC in the survey. 
In addition, Jordan et al. (2007; Paper XII) was able to 
determine the form of the GCLF for 89 galaxies in the 
survey. For galaxies where we were unable to measure 
the form of the GCLF we assume the mean to be the 
same as that of the cD galaxy VCC 1316, and a Gaus- 
sian distribution with sigma derived from Equation 2 in 
Jordan et al. (2006). In this paper, we will use the Gaus- 
sian parameterization of the GCLF. 

Most of our galaxies are small enough that the 
ACS/WFC field of view is sufficient to encompass all of 
the galaxy's GC system. In these cases, we determine the 
total number of GCs by counting GC candidates within 
our imaged field, correcting for incompleteness using the 
known mean and sigma of the GCLF. Calculating the 
completeness correction, however, is not straightforward 
because the surface brightness of the galaxy varies spa- 
tially. If, for instance, there were strong radial gradients 
in the GCLF or size distribution, we could be over- or 
under-correcting for unseen GCs. Our previous work on 
these galaxies (Jordan et al. 2005, 2007) shows that this is 
not likely to be the case. We use bright GCs (z < 22.5) — 
those that are bright enough as to be complete over the 
entire range of galaxy surface brightness and GC sizes — 
to sample the "true" distribution of background surface 
brightnesses and GC sizes. We then make the assump- 
tion that the fainter GCs sample the same distributions 
to derive a mean completeness correction for the entire 
population. 

In the more luminous galaxies, the ACS/WFC field of 
view is insufficient to image the entire GC system. To 
correct for this, we determine the radial spatial density 
profile of the GC system, fit a Sersic profile, and inte- 
grate over all radii to estimate the total number of GCs. 
For the largest galaxies, we supplemented the density 
profiles at larger radii with data from three sources: 1) 

pipctinc.html 
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Fig. 1. — Surface density of globular clusters in VCC 1316 
(M87/N4486) as a function of projected galactocentric radius. We 
use data from three sources to fit the spatial density profile: ACS 
pointings of the central field and companion fields (black dots), 
WFPC2 parallel pointings in adjacent blank fields (blue aster- 
isks), and ground-based data from McLaughlin (1999a) (orange 
diamonds). All data are in good agreement and are well fit by a 
Sersic profile (black line). We integrate the Sersic profile to obtain 
the total number of GCs. 

ACSVCS imaging of companion galaxies whose GC sys- 
tems are dominated by the halo population of the neigh- 
boring giant, 2) Our WFPC2 parallel imaging of nearby 
halo fields (Takamiya et al., in prep), and 3) Ground- 
based data from the literature. In cases where there is no 
wide-field ground-based data, our additional pencil-beam 
HST/WFC and WFPC2 observations still do a compara- 
bly good job because the spatial resolution of HST mini- 
mizes the noise from background contamination. We use 
the same method described above for completeness cor- 
rection, except we do so in bins of galactocentric radius 
instead of for the whole field of view. 

Errors for Nqc were determined using a Monte Carlo 
technique where we introduced the appropriate random 
Poisson noise to our density profiles, then fit and inte- 
grated the Sersic (1968) profile. After doing this 1000 
times, we used the distribution of estimated Nqc to de- 
termine the intervals that contained the closest 68% of 
the measurements to the mean. In Figure 1, we show the 
combined radial profile data for VCC 1316 (M87/N4486, 
hereafter M87), which includes the central ACS point- 
ing, four ACS pointings targeted at nearby compan- 
ions, WFPC2 pointings, and the ground-based data of 
McLaughlin (1999a). All are in good agreement and are 
well fit by a Sersic profile, except perhaps in the very 
central regions which are difficult to measure (especially 
from the ground) and do not contribute large numbers 
of GCs. The spatial density profiles of the GC systems 
of ACSVCS galaxies will be presented in Peng et al. (in 
prep), which will contain a more detailed description of 
the techniques. Below, we outline our use of supplemen- 
tary data in the larger galaxies. 

VCC 1226 (M49/N4472, hereafter M49). For the most 
luminous galaxy in the Virgo cluster, we used the GC 



counts in the fields of the neighboring compact galaxies 
VCC 1192 and 1199, excluding a region of 4 R e around 
the galaxies themselves, where R e is the effective radius. 
We also used four WFPC2 parallel fields in the halo, and 
the ground-based density profiles of McLaughlin (1999a) 
and Rhode & Zepf (2001). All these data are in excellent 
agreement. 

VCC 1316 (M87/N4486). Four ACSVCS galaxies are 
near enough to this cluster cD that the GCs observed 
in their images are dominated by those of the giant: 
VCC 1327 (NGC 4486A), VCC 1297 (NGC 4486B), 
VCC 1279 and VCC 1250. VCC 1250 also appears to 
have ongoing star formation. As described above for 
M49, we masked a region of 4 R e around these nearby 
neighbors and used the remaining GCs to constrain the 
density profile of the larger galaxy. In addition, we used 
five WFPC2 parallel fields, and the ground-based density 
profiles from McLaughlin (1999a) (see Figure 1). 
VCC 881 (N4406). We supplemented the profile with 
two WFPC2 fields and the ground-based density profile 
of Rhode & Zepf (2004). 

VCC 763 (N4374). We supplemented the profile with 
two WFPC2 fields and the ground-based density profile 
of Gomez & Richtler (2004). 

VCC 798, 731, 1535, 1903, 1632, 1231, 2095, 
1154, 1062, 1030, and 1664. For these galaxies, we 
extended the measured density profile using between 1 
and 3 WFPC2 fields per galaxy, which were at galacto- 
centric radii between 5'and 17'. 

When we compare the total number of GCs derived 
using the integrated GC radial density profiles against 
the total number counted within the ACS/WFC field of 
view, the numbers converge for galaxies with Mb > —18. 
Because integrated radial density profiles are more un- 
certain and not necessary in low luminosity galaxies, 
we have adopted Mb = —18 as a cutoff; for galaxies 
brighter than this, we use the Sersic-integrated GC num- 
ber counts, and for galaxies fainter than this, we use 
the corrected number of GCs directly counted within the 
ACS/WFC image. 

VCC 1938 has a close companion dE, VCC 1941, whose 
GC system complicates an accurate count of the larger 
galaxy's GC system. We mask a R = 70" region around 
VCC 1941, and count only GCs in the ACS/WFC field of 
view, so our count of the VCC 1938 GC system is likely 
to be a lower limit. 

All GC counts are corrected for foreground and back- 
ground contamination using control fields (except in the 
eight cases described below). In many other studies, the 
background is taken from an annulus around the target 
galaxy. The advantages of using separate control fields 
are: 1) We can sample a much larger area of sky to de- 
termine the mean surface density of contaminants, thus 
greatly reducing the Poisson errors introduced in contam- 
inant subtraction, 2) We can use the full field of view for 
the target galaxy to measure as much of the GC system 
as we can, 3) We can still study systems where the GC 
system fills the entire field of view, 4) We can naturally 
incorporate spatially varying completeness as a function 
of position within the galaxy. The only potential disad- 
vantage is the addition of cosmic variance into the error. 

The advantage in using multiple control fields over a 
local background approach is most clear for the dwarf 
galaxies that have few GCs, where we need to minimize 
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the error in the mean expected background. We deter- 
mine the error in the background by measuring the num- 
ber of GC-like objects selected in each of the 16 custom 
control fields, and taking the standard deviation of these 
counts. This number combines both Poisson errors and 
cosmic variance, and takes into account the varying se- 
lection and completeness from galaxy-to-galaxy. We find 
that the contribution of cosmic variance to the errors is 
on the order of the Poisson noise or less, and that the 
use of control fields is superior to the use of a local back- 
ground. 

For eight galaxies— VCC 1327, 1297, 1279, 1250, 1185, 
1192, 1199, 1178 — control fields are insufficient, and we 
measure the total number of GCs by counting candidates 
in a R = 70" aperture and subtracting a local back- 
ground. The first five are nearby neighbors to the giant 
elliptical M87 and so some (and sometimes all) GCs de- 
tected in their vicinity belong to the halos of the giants. 
The next three are neighbors of M49. The size of this 
aperture was chosen to fit the ACS/WFC field of view, 
and extends 11, 30, 4, 4, 10, 20, 11, and 5 R e in radius 
for these galaxies, respectively. These numbers are, in 
principle, lower limits to the total numbers of GCs in 
each system. 

In total, we present a homogeneous catalog of GC num- 
ber counts for 100 ACSVCS galaxies. We have taken 
great care in this paper to produce our best estimate for 
the total number of GCs in the ACSVCS galaxies. Peng 
et al. (2006a) listed numbers of GCs used in their analysis 
of color distributions, but we emphasize that those num- 
bers were not corrected for field of view or completeness 
and are best used only to evaluate the signal-to-noise ra- 
tio of the color distributions. Forbes (2005) and Miller 
& Lotz (2007) both derive Sn values from Peng et al. 
(2006a) (the Forbes (2005) paper uses them exclusively), 
and so the appropriate warnings apply. 

3.2. Total Luminosities and Masses of GCs 

Another quantity of interest is the total luminosity (or 
stellar mass) in GCs. We take a straightforward em- 
pirical approach to measuring the total luminosity. For 
each galaxy, we add up the z-band luminosity in ob- 
served GCs down to 1 mag fainter than the mean of the 
GCLF, with adjustments for contaminants and complete- 
ness. For galaxies with Mb < —18, we apply an aperture 
correction to the total luminosity in GCs, which is the 
ratio of the total number of GCs determined from the 
integrated radial profile to the total number observed in 
the ACS/WFC image. At this depth, we are directly 
counting 84% of the GCs in a Gaussian GCLF, but are 
sampling 99% of the luminosity in GCs. Thus, complete- 
ness corrections are small, and it has been advocated (e.g. 
Harris 1991) that the luminosity in GCs is a more robust 
quantity than the number of GCs. To obtain the total 
mass in GCs, we assume that each GC is a 13 Gyr simple 
stellar population (SSP) (Chaboyer et al. 1996) and use 
its (g-z) color in conjunction with the Bruzual & Chariot 
(2003; BC03) models with a Chabrier (2003) initial mass 
function (IMF) to obtain a mass-to-light ratio in the z 
bandpass. The range of M/L z variation with [Fe/H] for 
globular clusters is no more than ±20%. 

3.3. Galaxy Magnitudes: 
A Consistent Catalog of My And M z 



Globular cluster formation efficiencies can be measured 
against the total luminosity or mass of a galaxy. Spe- 
cific frequency is typically calculated using the absolute 

V magnitude. The Virgo Cluster Catalog lists Mb for 
all of our sample galaxies, but these magnitudes were 
not derived from CCD photometry. We do not have V 
imaging for our galaxies, so we predict V using optical 
colors measured from SDSS imaging. We use (ugriz) s d ss 
photometry measured directly from images in the Sloan 
Digital Sky Survey Data Release 5 (Adclman-McCarthy 
et al. 2007). In many cases the large sizes of the galaxies 
compared to the size of the SDSS CCDs required careful 
stitching of neighboring runs and camera columns to pro- 
duce flat images with matching sky (West et al. 2007). 
Using these specially prepared images, we measured total 
magnitudes in the five SDSS bands for all the ACSVCS 
galaxies using a growth curve analysis. Details of this 
analysis and the full catalog of magnitudes and colors 
will appear in Chen et al., in prep. 

For the purposes of this paper, we only use the SDSS 
photometry for two objectives, deriving a V magnitude, 
and supplementing the total z magnitudes given in Paper 

VI derived from ACS imaging. For determining each 
galaxy's V magnitude, we fit SSPs from BC03 to the 
four optical colors— (u - g) s dss, (g - r) sdss , (r - i) s dss, 
(i — z) s dss — and use the g s dss — V color of the best fit 
model. Across the entire sample, the g sds s — V color 
ranges from 0.18 to 0.48 with a mean of 0.37. In practice, 
this approach is extremely robust because g s d ss — V is 
almost entirely a function of (g — r) s d ss - We subtract this 
color from the measured g s dss magnitude to obtain V. 
The mean relationship between V and (g — r) s( [ ss matches 
the empirical relation of Blanton & Roweis (2007), except 
that our magnitudes are fainter in the mean by 0.03 mag. 

Because redder wavelengths are a better tracer of stel- 
lar mass than the traditionally used B and V, we will 
often use the total integrated z luminosities of galaxies in 
this paper. For most of these galaxies, we use the z mag- 
nitudes measured from the ACS images as presented in 
Paper VI. However, for the brightest 10 galaxies, the ACS 
field of view was substantially smaller than the extent of 
the galaxy, and we prefer to use the wide-field SDSS pho- 
tometry (after introducing a small (< 0.04 mag) color- 
dependent correction between the SDSS to the ACS pho- 
tometric systems). We also use the SDSS photometry 
for VCC 1030, 575, and 1512, which have suspect in- 
tegrated magnitudes from their ACS surface brightness 
profiles (see notes on these galaxies in Paper VI). Oth- 
erwise, the two independent measures of the total lumi- 
nosity are in good agreement. We choose to use the ACS 
magnitudes for the remaining 87 galaxies because they 
have higher signal-to-noise than the SDSS photometry, 
which were taken from the ground with much higher sky 
backgrounds and shorter exposures. 

3.4. Stellar Mass 

Stellar mass is a better basis for comparison when 
studying galaxies of different morphological types and 
star formation histories. The V band light in star- 
forming galaxies is affected significantly by young stars 
and is not the best tracer of the total mass. Using the 
(g-z) colors from Paper VI, J — K s colors from the Two 
Micron All Sky Survey (2MASS; Skrutskie et al. 2006), 
and the model SSPs from BC03 with a Chabrier (2003) 
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Fig. 2.— S N versus galaxy My for 100 ACSVCS galaxies. M87, 
the giant elliptical with Sn ~ 13 is a well-known outlier. Luminous 
early-type galaxies have higher Sn than intermediate-luminosity 
early-type galaxies by a factor ~ 2—3. Galaxies with intermedi- 
ate luminosities (—20.5 < My < —18) generally have Sn ~ 1.5. 
Early- type dwarf galaxies have a large spread in Sn, with some 
having zero GCs, and others having among the highest measured 
Sjv in our sample. The dotted line shows the mean trend (including 
M87) whose values are listed in Table 3. 

initial mass function, we obtain mean mass-to-light ra- 
tios in the z bandpass (A4/L z ) from which we can de- 
rive a total stellar masses. In the near-infrared, we mea- 
sure the J and K s galaxy magnitudes in images taken by 
2MASS, and supplement them with those in the 2MASS 
Extended Source Catalog (XSC) and Large Galaxy At- 
las (Jarrett et al. 2003). For 9 galaxies with unreliable 
K s photometry from 2MASS (those not included in the 
XSC), we use only the (g-z) color assuming a 10 Gyr 
SSP to determine M/L z . We assume a younger age than 
for GCs because dEs are measured to have younger ages 
than GCs or massive ellipticals (Geha et al. 2003). Al- 
though M./L is related to both age and metallicity, it 
is easier to determine A4/L than age or metallicity in- 
dividually. Most of the sample have M/L z consistent 
with old stellar populations (age > 5 Gyr), though a 
few galaxies are dE/dl transition objects and have no- 
ticeably younger mean ages and correspondingly lower 
A4/L z (the most extreme example being VCC 1499). 
For one of these galaxies, VCC 1030, we derive a low 
mean age, high metallicity, and a correspondingly low 
A4/L z . While this galaxy may be interacting, the col- 
ors may also be suspect due to a large-scale central dust 
disk. Until we get more data, we consider M/L for this 
galaxy uncertain. 

Fortunately, the z bandpass is less sensitive to recent 
star formation than bluer bandpasses, and the range of 
A4/L z is only 0.4 to 2.1, a dynamic range which is a 
factor of ~ 3 smaller than that for JvXJLb- The error 
in M/L z is dominated by the quality of the infrared 
photometry. We estimate errors in M/L z with a Monte 
Carlo procedure where the J — K s color is perturbed by 
a random amount drawn from a Gaussian distribution 
with sigma equal to the claimed photometric error. The 
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x Not including VCC 1316 (M87) 

A4/L z is then recalculated in this fashion 1000 times 
and the error is the half width of the middle 68% of 
the distribution. The typical error in Jvi/L z is ±0.4, or 
~ 25%. These errors are propagated forward into mass- 
normalized quantities, although any possible variations 
in the IMF are not included. 

All of the global galaxy properties discussed above are 
listed in Table 1, and all the quantities related to the 
globular cluster systems are presented in Table 2. 

4. RESULTS 

4.1. Sn versus Galaxy Magnitude 

The historical definition of specific frequency, Sn, is 
the number of GCs per unit My = — 15 of galaxy lumi- 
nosity, or 

S N = N GC X l 0-4(Mv + 15) (1) 

(Harris & van den Bergh 1981). This number is approxi- 
mately unity for the Milky Way. One of the motivations 
for calculating Sn is to understand whether GC forma- 
tion scales with the bulk of star formation in the same 
way across galaxy types and masses. In Figure 2 we show 
the behavior of Sn as a function of My in our sample of 
100 ACSVCS galaxies. The luminous early-type galaxies, 
mostly giant ellipticals, are well known to have Sn ~ 2- 
5. The cD galaxy, M87, has S N = 12.6 ± 0.8 which is 
consistent with other measures of its Sn (e.g., Harris 
et al. 1998). Early-type galaxies of intermediate lumi- 
nosity (— 22 < My < —18) have a nearly uniformly low 
(Sn) ~ 1-5. Galaxies in this luminosity range also have 
a tendency to be lenticular, with the VCC classifying 
70% of these ACSVCS intermediate-luminosity galaxies 
as SO, E/S0 or S0/E. By contrast, the fainter galaxies 
(My > —18) exhibit a wide range of Sn, with values as 
low as zero and as high as those of M87 19 . 



This trend has been hinted at previously from stud- 
ies of giant ellipticals, lenticulars, and dwarf ellipticals. 
In Figure 3, we show the ACSVCS sample combined 
with data compiled from the literature. At the high lu- 
minosity end, we take Sn values from the compilation 
of Ashman & Zepf (1998), using ones that were reli- 
ably determined from CCD data (gray diamond points). 

19 Strader et al. (2006) analyzed a subset of these data and 
claimed a bimodal distribution of Sn in the dEs. A histogram of 
the Sn distribution shows that it is better described as strongly 
peaked with a tail to higher Sn, especially when one includes dSOs. 
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Fig. 3. — 5jv versus galaxy My for 100 ACSVCS galaxies (black 
circles), and early- type galaxy data from the literature (references 
in the text). Literature values for Sn follow and extend trends 
visible in the ACSVCS galaxies. In particular, dwarf galaxy Sjy 
values from the work of Durrell et al. (1996) and Lotz et al. (2004) 
show that fainter dwarf ellipticals can have an even larger range 
of specific frequency. For clarity, error bars are not plotted for the 
HST/WFPC2 dEs, but the uncertainty can be very large. The 
dashed line shows the Sn value if a galaxy at that magnitude had 
1 GC. 

These include estimates from Kissler-Patig et al. (1996, 
1997), Dirsch et al. (2003a, 2003b), Dirsch, Schuberth 
& Richtler (2005), Forbes et al. (1996), Rhode & Zepf 
(2004), Zepf et al. (1995), and Harris, Harris & Geisler 
(2004). At the faint end, we include data from three 
sources which extend the range in galaxy luminosity by 
many magnitudes. We include the 7 dEs studied by Dur- 
rell et al. (1996a, 1996b) that do not overlap with our 
sample, the Virgo and Fornax dEs in the HST/WFPC2 
study of Miller et al. (1998) and Lotz et al. (2004), and 
the hve Local Group dwarfs listed in Lotz et al. 

At the high luminosity end, the literature values have 
large scatter, possibly due to their heterogeneous nature. 
One of the benefits of the ACSVCS is that relative dis- 
tances between Virgo cluster galaxies are both small and 
measured in a homogeneous fashion. At low luminosi- 
ties, the studies are mainly in Virgo and Fornax so rel- 
ative distances are less of a problem. At luminosities 
fainter than the ACSVCS, the trend of higher Sn cou- 
pled with a larger range of Sn continues to the limit 
of the data. The ACSVCS sample not only has smaller 
errors, it fills the important luminosity regime between 
giants and dwarfs. Kundu & Whitmore (2001a, b) stud- 
ied galaxies in the magnitude range —23 < Mb < —16, 
similar to the ACSVCS, and found that ellipticals gener- 
ally had higher Sn than lenticulars, but they were only 
able to measure "local" Sn within the WFPC2 field of 
view. 

Ultimately, studies of GC specific frequency are about 
how GC formation scales with galaxy mass, using My 
as a proxy. However, a redder bandpass more faithfully 
traces stellar mass, and in Figure 4, we introduce Sn, z , 
which is defined identically to Sn except that the num- 



FlG. 4. — Sn,z versus galaxy M z . This figure shows similar 
quantities to those in Figure 2 except the number of GCs is nor- 
malized to the galaxy absolute magnitude in the z-band, with M z 
is plotted along the x-axis. This far-red bandpass is a better tracer 
of the total stellar mass than the traditionally used V-band. The 
dotted line shows the mean trend with values in Table 4. 

ber of GCs is normalized to an absolute z magnitude of 
M z = —15, and plotted against M z . Sn.z values are 
1.5 to 3 times smaller than S n because early- type galax- 
ies are red (and thus more luminous in z). The specific 
frequencies of giant ellipticals are adjusted more than 
those of the dwarfs because massive galaxies are redder 
in color. While the trends in this figure are the same as 
those in Figure 2, one interesting difference is that the 
highest Sn.z values for the dwarfs now equal or exceed 
that of M87, with the change being due to M87's redder 
color. For the rest of the paper, we will use quantities 
based on M z . 

4.2. Sn,z and Bulge Luminosity 

In most scenarios of GC system formation, the GCs 
are associated with the stellar spheroids — either halo or 
bulge — and not with the formation of present day disks. 
If this is the case, then perhaps normalizing GC numbers 
to total spheroid luminosity would be more fundamental 
than normalizing to total luminosity. In our ACSVCS 
sample, over half of our galaxies are morphologically clas- 
sified in the VCC as either E/S0, S0/E, SO, or dSO. Most 
of these galaxies are at intermediate luminosities, and 
make up a substantial fraction of the galaxies that have 
low Sn, z - Could it be the case that these low Sn,z val- 
ues are due to the inclusion of a stellar disk in the total 
luminosity? Moreover, the scatter in the Sn.z of the in- 
termediate luminosity galaxies in our sample exceeds the 
scatter expected from the errors, implying that a param- 
eter such as bulge fraction might be important. 

We can test whether normalizing by bulge luminosity 
is more fundamental with quantitative bulge-disk decom- 
positions (Chen et al., in prep). In this paper, however, 
we use the morphological T-types given in the Third Ref- 
erence Catalog of Bright Galaxies (RC3) (de Vaucouleurs 
et al. 1991). These T-types are correlated with the ratio 



Formation Efficiencies of Globular Clusters 



9 




Fig. 5. — Sn,z versus galaxy M z for ACSVCS galaxies with VCC 
morphological types of E/SO, SO/E, SO, and dSO, and which have 
morphological T-types in the RC3. We plot the total Sjv % (same 
as Figure 4) as well as an arrow connecting to the galaxy's Spf z 
normalized to bulge luminosity alone. Some galaxies have SV,z 
similar to those in luminous ellipticals, but others are relatively 
unchanged. 

of bulge to total luminosity using the relation of Simicn 
& de Vaucouleurs (1986). We can do this for 55 galax- 
ies in the ACSVCS sample and we show how their «Sjv, z 
values change in Figure 5. This figure shows Sn, z as 
derived using total luminosity compared to Sn,z as de- 
rived from bulge luminosity alone with the two values 
connected by an arrow. For some galaxies, their specific 
frequencies do increase to values typical of the luminous 
ellipticals (SV,z~ 1 _ 2), but for others whose disks are 
less prominent, their Sn jZ stay relatively unchanged. For 
two galaxies, their RC3 morphological types values imply 
that they have little or no bulge, which give them very 
high S N ,z- 

Based on the current morphological classifications, it 
is not apparent that GCs form at a constant efficiency 
with respect to bulge luminosity. Normalizing to bulge 
luminosity increases the scatter in Sjv.z 111 the intermedi- 
ate luminosity range. However, a more thorough analysis 
based on modern bulge-to-disk decompositions is neces- 
sary before we can reach any strong conclusions. 



4.3. Normalizing to Stellar Mass 

Using the stellar masses for the ACSVCS galaxies, we 
calculate the T parameter introduced by Zepf & Ashman 
(1993), which is the number of GCs {N GG ) per 10 9 A4 Q , 

T = N GC /(M G */1O 9 M ) (2) 

where M G * is the stellar mass of the galaxy. The advan- 
tage of using T instead of Sn is that it allows compar- 
isons across galaxies with different mass-to-light ratios. 
In Figure 6a and b, we show T plotted against M z and 
Mg*- Although the errors are larger than in Figure 4, 
the same trends are evident. Our values for T are higher 



than those in previous studies, such as Rhode & Zepf 
(2004). This is due to differences in the mass-to- light 
ratios used. Previous studies have assumed M/Ly = 10 
for elliptical galaxies, whereas we have estimated M/L 
using galaxy colors and the BC03 models with a Chabrier 
(2003) initial mass function (IMF). Even accounting for 
the different bandpasses used, our M/L are systemati- 
cally lower. For example, for M49, the most luminous el- 
liptical in the sample, (V-z) = 0.97 and its M/L z = 2.0 
translates to M/Ly = 5, which is a factor of two lower 
than the canonical Zepf & Ashman values. We feel that 
although the M/L values in Zepf & Ashman (1993) have 
been valuable as standard conversions, they are too high, 
and that our values of M/L are more reasonable given 
recent dynamical mass measurements of elliptical galax- 
ies (Kronawitter et al. 2000; Cappellari et al. 2006). The 
use of a Salpeter (1955) IMF increases M/L by a factor 
of 1.8, but is not observationally motivated for stars be- 
low 1 Mq and would produce stellar M/L values higher 
than the total M /L measured in some early- type galax- 
ies. For the early- type galaxies in our sample, M/Ly 
ranges from 1 to 5, although some galaxies with obvious 
star formation have M/Ly < 1. We reiterate, though, 
that a using redder bandpass is better, and all our masses 
are derived using M/L z , and we discuss the V bandpass 
only for comparing with previous work. 

For the purposes of this paper, however, the absolute 
scale of T is not as important as the relative scale within 
the sample. Figure 6b shows that T spans a wide range at 
masses Mg* < 4xl0 9 A / f©, stays constant until IO^Mq, 
and then increases again at higher mass. 

4.4. Specific Luminosity and Mass 

If we want to know what fraction of a galaxy's lumi- 
nosity (or mass) is in GCs, we can total it up directly, as 
described in Section 3.2. We use the definition of specific 
luminosity, Sl presented in Harris (1991) except that we 
use the z bandpass instead of V: 



S L ,z = 100 x L GC \z/L gal 



axy,z 



(3) 

where Lqc,z and L ga i axy ^ are the total z-band lumi- 
nosities of the GCs and the galaxy. This quantity has 
two advantages over Sn in that it is independent of dis- 
tance and that the completeness corrections for unob- 
served faint GCs are extremely small. We also plot the 
specific mass, Sm, defined as 



S M = 100 x M GC /M G * 



(4) 



where M GG is the total stellar mass in GCs, calculated 
as described in Section 3.2. 

Figure 6c shows the specific luminosity Sl of GCs plot- 
ted against M z , and Figure 6d shows Sm versus M G *. 
The dashed line marks the 0.26% "universal" GC forma- 
tion efficiency from McLaughlin (1999a), and is a rea- 
sonable description of the GC mass fraction in inter- 
mediate mass galaxies. For the ACSVCS galaxies with 
-22 < M z < -19, (S L ) = 0.20, a s , = 0.14. Similarly, 
for galaxies with M G * = 0.4-6 x 10 W M Q , (S M ) = 0.17, 
°~Sm = 0.11. That the mass fraction is slightly lower than 
the luminosity fraction can be explained by the bluer col- 
ors and hence lower M/L of the GCs as compared to 
their host galaxies. 
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TABLE 4 
Bins of M z 
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Fig. 6. — Tjsr versus M z (a, upper left), and versus galaxy stellar 
mass (b, upper right). Tjv is the number of GCs per 1O 9 A^0. 
In the bottom two plots we show the specific z luminosity, Sl, 
versus M z (c, lower left) and the specific mass, Sm> versus stellar 
mass (d, lower right). The dashed line in (d) marks Sm = 0.26, 
the "universal" GC formation efficiency from McLaughlin (1999a). 
Because the ACSVCS sample is all early-types, A4/L does not vary 
much, and hence in all four panels we see trends similar to those 
for Sn and Sn,z ( see t ex t for exceptions). The difference between 
dwarfs and giants in Sl and Sm is n °t as large as it is in Sn 
or Tjv- This reflects the changing GCLF across galaxy luminosity 
(see Figure 7). In all figures, the solid line shows the mean trend, 
values in Table 5. 

The major difference in Figure 6 between the top two 
panels showing T (a and b) and the bottom two show- 
ing Sl, z an d Sm (c and d) is the comparison between 
the giants and dwarfs. While the dwarfs have very high 
number fractions, their luminosity and mass fractions are 
substantially lower compared to the giants. Why should 
this be the case? The reason for this has to do with 
the GCLF. If the GCLF was constant across all galaxies 
(i.e., same mean \i and width a), then Sl would mirror 
T. However, Jordan et al. (2006, 2007) showed that the 
GCLF varies as a function of galaxy mass in the sense 
that less massive galaxies host GC systems with fainter 
\x and smaller a. This has the effect of lowering the mean 
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Fig. 7. — Mean z luminosity of GCs, (L G c,z), in ACSVCS 
galaxies versus M z . If the GCLF was constant across all galaxies, 
(Lqq z ) would also be constant. The affect of a fainter GCLF 
turnover in dwarf galaxies (dotted line) only partially explains the 
lower mean luminosities. A combination of fainter turnovers and 
narrower GCLFs (dashed line) reproduces the observed trend. 



TABLE 5 
Bins of Mg* 



M* range 1 



(0.1,0.5) 


0.3 


18.2 


1.12 


0.30 


(0.5,2.2) 


1.0 


8.5 


0.51 


0.39 


(2.2, 10.0) 


4.8 


4.1 


0.21 


0.25 


(10.0,46.4) 


21.1 


4.7 


0.18 


0.17 


(46.4,215.4) 


80.8 


5.8 


0.21 


0.29 


(215.4, 1000.0) 
(215.4, 1000.0) 2 


321.0 


11.4 


0.38 


1.12 


324.7 


14.0 


1.00 


0.80 



1 M*/10 9 M Q 

2 Not including VCC 1316 (M87) 
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Fig. 8. — Fraction of red GCs (/ re£ j) versus M z . The dotted 
line represents the mean trend (Table 4). More luminous galaxies 
have higher fractions of red GCs up to M z ~ —22, but at higher 
luminosities there is a flattening or turnover in / ret j, and the most 
luminous galaxies do not have the highest red GC fractions. 

luminosity (and mass) of GCs in the dwarf galaxies. We 
illustrate this in Figure 7, where we plot the mean z lu- 
minosity of GCs in each galaxy, (Lac, z ), against galaxy 
luminosity (in this case Mb because that is how we pre- 
sented the data in Jordan et al.). It is clear that a con- 
stant GCLF does not describe the data well. The dot- 
ted line shows the effect of varying /i, and the dashed 
line shows the relation predicted by the combination of 
fainter /i and a narrower a in the dwarfs. Although the 
scatter is large, the expected change in the GCLF can 
account for the changes in mean GC luminosity that we 
see. This trend is actually more pronounced in mass than 
in luminosity because GCs in dwarfs are more metal-poor 
and have lower M/L than those in giants. 

The result is that some of the more extreme Sn or 
T values seen in the dwarfs are somewhat less extreme 
when expressed as a mass fraction. Nevertheless, many 
low-luminosity galaxies still have Sl and Sm that are sig- 
nificantly higher than those in the intermediate- L galax- 
ies. Because the global trends are similar in Sn, T, Sl, 
and Sm, we will refer to these quantities collectively as 
"GC fractions" . 

4.5. Specific Frequencies of Red and Blue GCs 

Our previous study of GC color and metallicity distri- 
butions in the ACSVCS galaxies show that they are, on 
average, either bimodal or asymmetric across the entire 
luminosity range of the sample, and we use the prod- 
ucts of that analysis (Peng et al. 2006a) in the current 
study. The blue (metal-poor) and red (metal-rich) are 
believed to trace either different epochs of formation or 
different progenitor halos (however, see Yoon, Yi & Lee 
2006 and Cantiello & Blakeslee 2007 for arguments that 
the bimodality may be an observational consequence of a 
nonlinear metallicity-color relation). It is an interesting 
question to ask how the specific frequencies of each GC 
population scale with galaxy properties. For example, 



Rhode & Zepf (2004) showed that the mass-normalized 
number of blue GCs, Tu U e, increases as one goes from 
spirals to elliptical galaxies, although their spirals were 
less massive than ellipticals so the trend could also have 
been one in galaxy mass. Assuming that mergers of spi- 
rals only produce new red GCs, they argued that the GC 
systems of ellipticals cannot be formed purely by spiral 
mergers. 

We determine the fraction of blue and red GCs in a 
hybrid approach similar to what we do for total num- 
bers. For the brightest 21 galaxies in the sample (a 
complete magnitude-limited sample, as ranked by Bt in 
Cote et al. 2004; Paper I) there are sufficient numbers 
of GCs that the KMM two-Gaussian fits to the (g-z) 
distribution performed by Peng et al. (2006a; Paper IX) 
are reliable enough that we can use the "dip" between 
the two Gaussians — the color at which a GC is equally 
likely to belong to the blue or red GC distribution — as 
the dividing color between the two populations. For the 
remaining galaxies, we assume a fixed dividing color of 
(g-z)= 1.16. Although the colors of the individual peaks 
vary as a function of galaxy luminosity, the dip color is 
relatively invariant (see Figure 5 in Peng et al. 2006a). 

The brighter galaxies may have better defined color 
distributions, but we only observe a fraction of their en- 
tire GC system, and the red-to-blue ratio is observed to 
decrease as a function of galactoccntric radius. We cor- 
rect for this bias in the brightest 14 galaxies by fitting 
Sersic models to the surface density profiles of the red 
GCs separately, and integrating the best model to ob- 
tain their total numbers. We fit to profiles derived from 
our ACS and WFPC2 data, as described in §3.1. In the 
most luminous galaxies such as M49 and M87, the frac- 
tion of red GCs within the ACS/WFC is as high as 0.6, 
but the red fraction of the entire GC system is more like 
0.3. For the remainder of the galaxies, we find that the 
ACS/WFC encompasses nearly all of the red GC system 
and use our corrected counts to determine the red GC 
fraction. 

Figure 8 shows the fraction of red GCs as a function 
of M z . The fraction of red GCs, j re d, generally increases 
with galaxy luminosity, going from ~ 0.1 to ~ 0.5, but 
at M z < —22, the trend in f rec i appears to either flatten 
or perhaps turn over. This is in contrast to the results 
of Paper IX which only quantified the red GC fraction 
within the ACS/WFC images, and thus were biased to 
detecting the more centrally concentrated red GCs. With 
the proper aperture corrections, we can see that the most 
luminous galaxies in the sample do not have increasingly 
higher fractions of red GCs. Is this because there are 
fewer red GCs in these galaxies or more blue GCs? 

Figure 9 shows the Sn, z for blue and red GCs as func- 
tion of M z . The trend in the blue GCs mirrors the overall 
trend seen in Figure 4, with the massive and dwarf galax- 
ies having the highest Sn, z - This is not surprising on the 
faint end since most of the GCs in the fainter galaxies 
are blue. Even for the massive galaxies, the high Sn, z 
values are dominated by blue GCs. However, the spe- 
cific frequencies of red GCs also exhibit an increase for 
M z < —21, especially for the cD galaxy M87, which is 
a 4cr outlier. The elevated Sn, z ,Muc and SN, z ,red for the 
most massive galaxies, suggests that massive galaxies are 
not underproducing red GCs, which is one possible in- 
terpretation of Figure 8. Instead, these massive galaxies 
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Fig. 9. — Sjv.z versus M z for blue (left) and red (right) globular 
clusters. Both red and blue GCs show enhanced Sjv, z in massive 
galaxies, but the variation in Sjv z across galaxy mass is dominated 
by the blue GCs. This is true even in the massive galaxies. The 
dashed lines show the mean trends with values listed in Table 4. 

have more red GCs, but even more blue GCs. 

4.6. Sn.z and Nucleation 

Observations by the HST/WFPC2 of dEs in the Virgo 
and Fornax Clusters (Miller et al. 1998; Lotz et al. 2004; 
Miller & Lotz 2007) have found that dEs with stellar 
nuclei have a higher mean specific frequency, possibly 
implying that a higher past star formation efficiency re- 
sulted in the formation of both nuclei and GCs. The dEs 
in our sample are more luminous than the ones studied 
in the WFPC2 snapshot survey, and almost all of the 
galaxies are nucleated (Cote et al. 2006; Paper VIII). 
Therefore, it is difficult to test if there is any correlation 
between Sn,z and nucleation. Only four dwarf galaxies 
are definitely non-nucleated (Type II in Table 1 of Pa- 
per VIII)— VCC 1049, 1833, 1499, and 1512. Of these 
galaxies, VCC 1499 and 1512 are dE/dl transition ob- 
jects with young stars, and VCC 1049 has bluer colors 
toward the center. The Sn z values of these galaxies 
range from 0.97 (VCC 1049) to 3.14 (VCC 1499) and do 
not appear different from the rest of the sample. When 
normalized to stellar mass, however, the young stellar 
populations have lower A4/L and thus much higher T 
and Sm values. 



4.7. Sn,z and Environment 

One of the more intriguing questions presented by 
plots such as Figure 4 is the nature of Sn,z in the low- 
luminosity galaxies (M z > —19), which we will loosely 
refer to as "dwarfs" , although some of the galaxies in our 
sample are on the more massive end of the spectrum of 
dwarf galaxies. Our relatively small errors show that the 
large spread in specific frequency for dwarf galaxies is not 
simply due to the expected observational scatter. Given 
our errors, the observed distribution of Sn.z is 2.4 times 
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Fig. 10. — Sjv.z vs projected clustercentric distance for low- 
luminosity galaxies (M z > —19), excluding VCC 571 and 538 
which are known to be ~ 6—7 Mpc behind the Virgo core (Mei 
et al. 2007). The center of the Virgo cluster is taken to be the lo- 
cation of the cD galaxy, M87, and the top axis shows R P /r200 where 
r 200 = 1-55 Mpc. There is a notable trend in specific frequency 
with clustercentric radius. All but one galaxy with Sjv.z > 2 are 



TABLE 6 



R p range 




Sn,z 


T 


£>N,z,blue 




Mpc 


Mpc 




(0.00,0.15) 


0.07 


0.47 


3.9 


0.40 


0.07 


(0.15,0.30) 


0.23 


3.46 


33.3 


3.17 


0.29 


(0.30,0.50) 


0.40 


1.58 


15.7 


1.34 


0.23 


(0.50, 1.00) 


0.83 


1.38 


13.0 


1.14 


0.23 


(1.00, 1.50) 


1.24 


0.70 


8.4 


0.59 


0.11 


(1.50,2.00) 


1.62 


1.05 


11.6 


0.99 


0.06 



broader than we would expect if all dwarfs had a single 
Sn,z ■ This implies that there is at least one other param- 
eter besides galaxy mass that governs the GC formation 
efficiency. 

In this section, we investigate the relationship be- 
tween specific frequency and the galaxy's environment — 
specifically, its distance from the cluster center, which is 
taken to be the location of the cD galaxy, M87 (Binggeli 
et al. 1985). Environment clearly plays a role in galaxy 
evolution, as is evidenced by the morphology-density re- 
lation, and this can be through gravitational and hydro- 
dynamic processes such as tidal and ram pressure strip- 
ping, or through initial conditions where halos in denser 
regions collapse earlier. Dwarf galaxies, being the most 
vulnerable, are more likely to express the effects of their 
environment. 

Figure 10 shows Sn.z for dwarf galaxies against their 
projected distance from the cluster center (R p ), ignoring 
VCC 571 and 538 which are known from their SBF dis- 
tances to lie 6-7 Mpc behind the cluster core (Mei et al. 
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Fig. 11. — 5jv versus three dimensional clustcrcentric distance for 
low-luminosity galaxies (M z > —19). The top axis shows R 3 d/ t 200 
where r200 = 1-55 Mpc. The error bar at top right shows the mean 
distance error. Points with arrows do not have SBF distances so 
we have used their projected distances and added 0.3 Mpc (the 
median (R 3 d — Rp) for the rest of the sample). The left error bar 
is at their projected radii, and are thus lower limits on their 3-d 
radii. 
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Fig. 12. — Spatial distributions of GC candidates (circles) in 4 
dEs with the highest specific frequencies (VCC 1407, 1545, 1539, 
and 230). In all cases, there is evidence of a rich GC system that 
is centrally clustered around the galaxy, showing that their ele- 
vated Sn is not due to an enhanced level of interloping GCs from 
the cD galaxy or an intracluster population. These galaxies have 
~ 5x more GCs than "normal" galaxies with Sjv, z = 1- The images 
show the entire ACS field of view, and the scale bar in the lower 
left of each image has a length of 30". The blank diagonal strip 
is the gap between the two ACS CCDs. Catalogs have been sta- 
tistically cleaned using expected contamination determined from 
control fields. 



TABLE 7 

Bins of R 3d for Galaxies with M z > -19 



R 3d range 


(Rsd) 


£>N,z 


T 


^N,z,blue 


^N,z,red 


Mpc 


Mpc 








(0.00,0.25) 


0.19 


1.38 


11.3 


1.19 


0.18 


(0.25,0.50) 


0.34 


2.39 


23.9 


2.19 


0.20 


(0.50, 1.00) 


0.71 


2.13 


20.8 


1.87 


0.26 


(1.00, 1.50) 


1.27 


1.27 


13.5 


1.05 


0.22 


(1.50,2.00) 


1.65 


0.78 


7.7 


0.68 


0.10 


(2.00,2.50) 


2.14 


0.50 


5.2 


0.47 


0.03 



2007). There is a remarkably clear correlation in which 
dwarfs closer to M87 have higher Sn, z - Of the 14 galaxies 
with Sn^ z > 2, 13 lie within a projected radius of 1 Mpc 
(the exception being VCC 21, a possible dE/dl transi- 
tion galaxy). We note that this is roughly half of r^m 
for the Virgo A subcluster centered on M87 (McLaugh- 
lin 1999b; Cote et al. 2001), where r 200 = 1-55 Mpc 
is the radius at which the mass density of M87 is 200 
times the critical density. This value has been recal- 
culated for Dm87 = 16-5 Mpc using the McLaughlin 
(1999a) model (D. McLaughlin, priv. comm.). When 
we plot Sn,z against three-dimensional clustercentric dis- 
tance (Figure 11), we see the same trend although the 
error in the line-of-sight distance introduces more scat- 
ter. Tables 6 and 7 list the values for the mean trends in 
these Figures. 

The M87 globular cluster system is one of the most 
extreme in the local supercluster and it is possible that 
some of the M87 halo GCs, or a population of intra- 
cluster GCs, are contaminating the GC systems of the 



nearby dwarfs. In fact, the five galaxies closest to the 
cD, VCC 1297, 1327, 1279, 1250, and 1185 have been 
treated differently in our analyses in this paper because 
it is obvious that some to all the GCs detected are 
part of the M87 halo. One of the telltale signatures of 
this is whether the GCs are uniformly distributed across 
the ACS/WFC image rather than centrally concentrated 
around the targeted galaxy. Figure 12 shows the loca- 
tions of GC candidates around the 4 dEs with the high- 
est 5jv, z - In all cases, the star clusters are concentrated 
toward the center of the galaxy and not in a uniform 
spatial distribution. An extension of the M87 GC radial 
density profile using a Sersic model fit to the ACSVCS 
data and the ground-based data of McLaughlin (1999) 
(see Figure 1) predicts only 1-3 M87 GCs over the entire 
ACS/WFC field of view at the distances of VCC 1407, 
1545, and 1539. The GC radial profile from Tamura et al. 
(2006) gives slightly higher values, with an expected 9^g 
M87 GCs in the ACS field at the projected distance of 
VCC 1407, a dE which has 50 GCs. Given the difficul- 
ties in measuring the M87 GC density from ground-based 
data, especially given uncertainties in background sub- 
traction, these independent estimates are consistent. Ex- 
trapolations of a de Vaucouleurs profile fit to the Tamura 
et al. data predict - 2 M87 GCs in the ACS fields for 
VCC 1545 and 1539. We also have a parallel WFPC2 
observation at a distance of 40' from M87 in which we 
find a number of GCs consistent with zero. Finally, as 
these galaxies have 31-54 GCs each, this contamination 
from M87 GCs is expected to be at most 6-18%. 

The general trend of dEs having higher Sn z 1S reversed 
for the two dwarfs (VCC 1297 and 1185) with R p < 
100 kpc. These galaxies have low intrinsic numbers of 
GCs and we hypothesize that they have had their GCs 
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Fig. 13. — Spatial distribution of ACSVCS galaxies color coded 
by S jy z . Specific frequency increases as colors change from blue to 
red. Filled points are dwarfs (M z > —19) and open points are the 
more luminous galaxies. The three large solid circles represent the 
three massive ellipticals in the Virgo cluster: the cD galaxy M87 
center, M49 bottom, and M60 (VCC1978/N4649) left. The large 
dotted circle represents r20o/2 = 775 kpc for the M87/ Virgo A 
subcluster. The high Sjv, z galaxies are preferentially around the 
cD, while no enhancement appears around the cluster's most lumi- 
nous galaxy M49. Dwarfs immediately in the vicinity of the giants 
have low Sjv z arid may be tidally stripped. 

stripped from them by M87. Both VCC 1297 and the 
more luminous VCC 1327 are within R p i=s 40 kpc, and 
have GC numbers consistent with zero. The galaxies 
VCC 1192 and 1199, which are at comparable distances 
from M49, also have undetectable intrinsic GC systems. 

Figure 13 shows the positions of the ACSVCS galaxies 
on the sky, color and size coded by Sn, z , with open points 
representing giants and filled points for dwarfs. The high- 
Sn,z dwarfs are clustered around M87, and there is no 
similar effect around the brightest galaxy in the cluster, 
M49, nor around VCC 1978 (M60/N4649). Interesting 
exceptions are the aforementioned galaxies within the im- 
mediate vicinities of the giants (R < 40 kpc), which have 
few or no GCs and may have been tidally stripped of their 
GC systems. The specific frequencies of the more mas- 
sive galaxies are generally uniform, as shown in Figure 4, 
and thus are relatively unaffected by their distance from 
M87. 

We can re-plot Figure 4, except divided into two sam- 
ples, one inside and one outside a projected clustercentric 
distance of 1 Mpc. Figure 14 shows the difference be- 
tween the two samples. Galaxies within 1 Mpc mirror 
the trends seen in Figure 4, with both giant and dwarf 
galaxies having high Sn.z- However, the sample of galax- 
ies outside of 1 Mpc entirely lacks high-SV^ dwarfs. The 
division is not perfect, as there are still \ow-Sn.z dwarfs 
within 1 Mpc (some of which have larger 3-D distances), 
but the lack of high-5Vz dwarfs on the outskirts of the 
cluster or around the other massive ellipticals is very 
clear. The one possible exception to the general trend 
is VCC 230, which lies at R p = 0.96 Mpc and has an 
Sn,z higher by a factor 2 than other galaxies at compa- 



FlG. 14. — Sm,z versus galaxy M z , as in Figure 4, except divided 
by clustercentric distance. Galaxies within 1 Mpc of the cluster 
center exhibit the full range of Sjvzi but nearly all the dwarfs 
outside of 1 Mpc have low specific frequencies (Sm,z < 1-5). The 
only exception is a dE/dl transition object (VCC 21). 




R p (Mpc) 

Fig. 15. — Sm z versus projected clustercentric distance for blue 
GCs (left) and red GCs (right) in low-luminosity galaxies (M z > 
— 19). Given that most GCs in dwarfs are metal-poor, it is not 
surprising that the trend for blue GCs mirrors that in Figure 10. 
The Sjv,z for the red GCs may also show a tendency to be enhanced 
at small clustercentric radii, but the numbers involved are small. 

rable distances. 

The GCs in dwarf galaxies are predominantly from 
the blue (metal-poor) subpopulation, but many in the 
ACSVCS sample also possess small numbers of red GCs. 
Are both of these subpopulations affected by environ- 
ment? In Figure 15, we plot Sjf x in dwarf galaxies 
against R p for the blue and red GCs separately. As ex- 
pected, the blue GCs mirror the overall trends since they 
dominate the GC budget. For the red GCs, although the 
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Fig. 16. — Color-magnitude properties of dwarfs with high and 
low Sjv z- We show (g—z) vs. M z and accompanying histograms 
for the ACSVCS sample of galaxies. Galaxies with M z > —19 
are separated at Sn z = 2. The luminosity distribution of the 
two subsets are nearly identical, while the high-Sjv.z dwarfs have 
a median color that is formally redder than that of the low-Sjv,z 
dwarfs (shown by dashed lines in right hand histogram). 

errors are large, there is a slight hint that dwarfs with 
smaller clustercentric distances tend to have elevated red 
GC specific frequencies. The mean SN, z ,red of dwarfs 
within 1 Mpc is 0.2, compared to 0.1 for those outside. 

Do the dwarfs with higher Sj? z exhibit any intrinsic 
characteristics that differentiate them from the others? 
Figure 16 shows the (g—z)-M z color-magnitude diagram 
for the ACSVCS galaxies, singling out the dwarfs with 
both low and high Sn.z, and dividing the sample at 
Sn,z= 2. The top histogram plots the M z distribution 
of the total, high-S^r 2 , and low-Sjv,z samples and shows 
how the luminosity distributions of the two groups are 
nearly identical. On the right-hand histogram, we plot 
the color distribution with the median colors marked by 
the dashed lines. Although the high-«Sjv,2 dwarfs formally 
have redder colors, A(g — z) = 0.05 mag, the difference is 
not significant. We have also investigated possible differ- 
ences in their structural properties (Sersic n, r e ) as well 
as in the colors of their globular clusters, and for none of 
these properties are high-SV^ dwarfs significantly differ- 
ent from low- Sat. z dwarfs. 

5. DISCUSSION 

5.1. Mass Dependence of the GC Mass Fraction 

Globular cluster specific frequency and its related 
quantities are clearly dependent on the mass of the host 
galaxy. However, GC fraction does not vary monoton- 
ically with galaxy mass, making it unlike most other 
properties measured for the ACSVCS galaxies — e.g., col- 
ors, structural parameters (Ferrarese et al. 2006a), core 
deficit and excess (Cote et al. 2007), GC mean colors 
(Peng et al. 2006a), GC mean sizes (Jordan et al. 2005), 
and GC luminosity function parameters (Jordan et al. 
2007). These papers (particularly Papers VI and VIII) 
show that there is no clear distinction between "dwarfs" 



and "giants" , and that early-type galaxy properties have 
a smooth, monotonic dependence on galaxy mass. For 
GC fractions, however, a transition does appear to exist. 

It is possible that GC fraction is telling us more about 
the variable formation efficiency of stars in the field than 
it is about the formation of GCs themselves. Blakeslee 
et al. (1997), Blakeslee (1999) and McLaughlin (1999a) 
have shown that, at least on the high mass end, the num- 
ber of GCs appears to scale directly to the baryonic, or 
even the total mass. Kravtsov and Gnedin (2005) show 
a similar relation in their simulations of GC formation. 
High Sn values may then be interpreted as a lower frac- 
tion of baryons that form field stars. 

The mismatch between total mass and stellar mass 
across the galaxy mass function is now well known in 
the study of galaxy formation. Observationally, dynam- 
ical mass estimates of galaxies across a wide range of 
mass require high mass-to- light ratios for both high (Cote 
et al. 2001, 2003) and low mass galaxies (Mateo 1998), 
with low Ai/L measured for galaxies around L* (Ro- 
manowsky et al. 2003; Peng et al. 2004; Napolitano et al. 
2005). Mass-to-light ratios for ensembles of galaxies de- 
rived from weak lensing have also produced similar trends 
in Ai/L (e.g., Guzik & Seljak 2002; Hoekstra et al. 
2005; Mandelbaum et al. 2006). One can also take a 
statistical approach and match the expected dark mat- 
ter halo mass distribution from simulations and the ob- 
served galaxy luminosity function (Berlind & Weinberg 
2002; van den Bosch et al. 2003, 2007; Vale & Ostriker 
2007), which does not require us to know the details 
of galaxy formation. These halo occupation studies in- 
fer a maximum conversion efficiency of baryons to stars 
at halo mass Aih ~ 2 x lO 11 ^©, or a stellar mass of 
M-G* ~ 7.5 x 1O 9 .M0. In simulations of galaxy forma- 
tion, this AAh/L-L relation (or alternatively, Aih/ M-g*~ 
A4g*) requires various heating mechanisms (photoioniza- 
tion, stellar and AGN feedback) to prevent the formation 
of more stars in galaxies than are observed (Benson et al. 
2002; Croton et al. 2006). 

What happens if we make the simple assumption that 
Nqc oc Mh, or equivalently Sm °c Aih/Aic*? Ideally, 
we would know the total dynamical mass of each galaxy 
in our sample. Without this information, we can apply 
average halo mass-to-light ratios derived from the halo 
occupation studies, in this case using the parameteriza- 
tion of van den Bosch et al. (2007) and a WMAP3 cos- 
mology (Spergel et al. 2007) . We transform their Lb to 
A4g* using a range of AAg*/Lb, monotonically increas- 
ing from 1.8 to 4.1 as a function of galaxy luminosity, 
derived from a polynomial fit to the inferred A4g*/Lb 
and measured L z of the ACSVCS galaxies. This results 
in a simple prediction for the behavior of Sm as function 
of galaxy stellar mass. 

In Figure 17, we plot the ACSVCS data for Sm against 
Aic* with errors for both Sm and A4g+- We extend the 
data to lower mass galaxies by including the Virgo and 
Fornax dEs and five Local Group dEs from Lotz et al. 
(2004). For these latter data, we assume A4g*/ Lb ~ 2.5 
(an extension of the mass-to-light ratios fitted to the 
ACSVCS dwarfs) and (M G c) = 1-6 x W 5 M Q , the 
mean GC mass in our ACSVCS dwarfs. The Virgo and 
Fornax dE points do not include errors in their stellar 
masses, and for the Local Group dEs, we have assumed 
an error in their luminosities of 0.25 mag. Overplotted 
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Fig. 17. — Mass fraction of GCs, >Sm> versus galaxy stellar 
mass M.Q* f° r ACSVCS galaxies (filled dots), Virgo and Fornax 
dEs from Lotz et al. (2004) (crosses), and Local Group early-type 
dwarfs (diamonds). We overplot the expected behavior of Sm 
assuming that the mass in GCs follows the total halo mass as in- 
ferred from the A4h/L-Mh relation of van den Bosch (2007). This 
assumption can explain the rise in Sm for luminous galaxies, al- 
though the cD galaxy M87 is off the relation because its stellar 
mass is not representative of its halo mass. At the low mass end, 
we also expect a sharp rise in Sm which is traced by some cluster 
and Local Group dEs. At masses of ~ 10 9 A4q, however, many dEs 
in the ACSVCS sample have higher Sm than might be expected. 
The high-Svn and \ow-Sm dwarfs cannot be explained simultane- 
ously if GC mass fraction is solely a function of galaxy stellar mass. 



on the data is the predicted behavior of Sm assuming 
Sm -M-h/ M-g*- The only truly free parameter in 
the plot is the vertical normalization, which is chosen 
to match both the massive elliptical M49 and the Lo- 
cal Group dwarf spheroidal galaxies Fornax and Sagit- 
tarius. The assumption that the GC mass fraction is 
proportional to the inverse of the stellar mass fraction, 
Sm k M-h/Mc* predicts that Sm should be high at 
both high and low masses. At high stellar masses, the rise 
reflects the increasing M/L of massive elliptical galax- 
ies. One galaxy, M87, is well off the predicted relation, 
but this is because its GC population is more likely rep- 
resentative of the entire Virgo A subclustcr within which 
M87 resides (McLaughlin 1999a). The sharp cutoff in 
the galaxy luminosity function and the nearly univer- 
sal luminosity of brightest cluster galaxies means that 
Sm can vary widely depending on whether the luminos- 
ity is associated with a galaxy halo or the cluster halo 
(Blakcslce 1999, Jordan et al. 2004a). At low masses, 
Sm is also predicted to rise steeply, but this is at masses 
lower than those in our sample. Many of the ACSVCS 
dwarfs around Mg* ~ lO 9 ^© have Sm higher than we 
would infer by a factor of 4-8, and these are the same 
ones that are in the inner regions of the Virgo cluster. At 
masses below Mg* ~ 10 8 Mq, the dEs from Lotz et al. 
(2004) are either consistent with or below the expected 
rise, although the observational error is quite large (and 
not shown for clarity). So, while the linear scaling of 
Mgc ex Mh does seem to be relevant to Sm across a 



large galaxy mass range, the Sm f° r the ACSVCS dwarf 
galaxies in this diagram is not explained by a dependence 
on galaxy mass alone. The scatter in Sm f° r dwarfs (and 
also possibly for intermediate mass galaxies) also can- 
not be explained by the intrinsic scatter in galaxy M/L 
(from van den Bosch et al. 2007). The scatter is even 
harder to explain for Sl where the errors are consider- 
ably smaller. 

Both Forbes (2005) and Bekki, Yahagi & Forbes (2006) 
have also attempted to explain the increase in Sn for 
dwarf galaxies by assuming a relation between GCs and 
galaxy mass. These attempts followed similar argu- 
ments by Durrell et al. (1996) and McLaughlin (1999a). 
Forbes (2005) used the data for ACSVCS galaxies pre- 
sented in Peng et al. (2006a, Paper IX) to compare 
the approximate Sn of blue GCs in dwarf galaxies with 
the scaling relation of Dekel & Birnboim (2006), where 
M/L v oc M~ a . Forbes (2005) use a = 2/3 from Dekel 
& Birnboim (2006) for masses below the critical galaxy 
stellar mass of 3 x lO 10 ^© where galactic winds should 
be more efficient at blowing gas and metals from the 
galaxy (Kauffmann et al. 2003; Tremonti et al. 2004). 
They also assume a constant number of GCs as a func- 
tion of galaxy mass, giving Sn oc Al~ 5 / 3 . Although this 
formulation predicts a rapid rise in Sn and Sm for fainter 
galaxies, it is much too steep to accommodate the data in 
the way that they present it. The normalization of this 
relation as presented in Forbes (2005) is not constrained, 
but because the relation is steep, the normalization is 
critical for comparing to the data. If we fix to the mean 
Sm at the critical mass, the relation overpredicts Sm f° r 
all the dwarfs. This is due to the fact that dwarf galax- 
ies do not have a constant number of GCs. Bekki et al. 
(2006) simulate GC systems and test various values for a, 
concluding that the observations are not consistent with 
a = (a constant M/L), and it would be interesting 
to determine a again with this new data set. Whatever 
the best value of a, however, a simple relation where the 
specific frequency is determined solely by the present day 
stellar mass of a galaxy is unable to describe all of the 
data on the low mass end. 

5.2. Environment and Biased GC Formation 
5.2.1. Cluster Dwarfs in the Millennium Simulation 

As we showed in Figures 10-15, GC specific frequency 
is not purely a function of galaxy mass, but can also 
depend on environment. We do not see environmental 
effects on Sn for massive galaxies with the exception 
of M87, whose location defines the cluster center. For 
intermediate-mass galaxies in particular, Sn does not 
appear to depend on environment. For dwarf galaxies, 
however, specific frequency is a strong function of prox- 
imity to the center of the Virgo Cluster. Nearly all dwarfs 
with high Sn are within R p = 1 Mpc. 

One of the main problems with building up the blue 
GC sub-population in the halos of massive galaxies is 
that the ratio of metal-poor GCs to metal-poor field stars 
in galaxies is very high (Harris & Harris 2002). In other 
words, if halos are built up through the accretion of dwarf 
galaxies or dwarf-mass fragments then their specific fre- 
quencies need to be high in order to keep the main body 
of the galaxy metal-rich. These central dwarfs with high 
Sn may then hold the key to the formation of GC sys- 
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Fig. 18. — The mass- weighted stellar age of low mass (M z > —19) 
early-type cluster galaxies in the Millennium Simulation. Galaxies 
in the inner regions of the galaxy cluster have older mean ages. The 
black error bars with hats depict the 2a error in the mean, while 
the gray error bars indicate the la width of the age distribution in 
each bin. 

terns, as they are likely to be most similar to the pro- 
genitors of any dwarfs that have since merged into the 
halo of M87. Figure 9 shows that the Sn.zMuc values for 
the innermost dwarfs are as high or higher than that of 
M87 or any of the other giant ellipticals. These dwarfs 
could be the survivors of an accreted population of pro- 
togalaxies that may have had even higher GC mass frac- 
tions and redder GCs. In fact, the two galaxies closest to 
M87— VCC 1327 and 1297— have GC systems that are 
either entirely stripped or undetectable against the M87 
GC system, and their g — z colors are much redder than 
those of other galaxies with the same luminosity (see the 
red outliers in Figure 16). 

What is the cause of the enhanced GC fractions in 
these central dwarfs? Are they more efficient at produc- 
ing GCs for their mass, less efficient at forming field stars, 
or better at keeping the GCs they form? As discussed 
earlier, one possible way for GC production to be biased 
towards the central dense regions of the galaxy cluster 
is for these galaxies to form a larger fraction of their 
stars earlier and at a higher star formation rate density. 
Although we do not currently have the ability to deter- 
mine the detailed star formation histories of the ACSVCS 
galaxies at these early times, we can use simulations to 
study global trends in the star formation histories of low 
mass early- type galaxies in the cluster environment. By 
using simulations, we can test the consistency of the hy- 
pothesis that the central dwarfs must be both older and 
have higher SFR densities. 

The simulation we use for comparison is the "Millen- 
nium Simulation" carried out by the Virgo Consortium 
(Springel et al. 2005) coupled with the semi-analytic 
models presented by De Lucia et al. (2006). This simu- 
lation consisted of 2160 3 dark matter particles followed 
from z = 127 to the present day in a volume 500/i _1 Mpc 
on a side. The spatial resolution is 5/i _1 Mpc, and the 



simulation is essentially complete for all galaxies with 
stellar mass greater than 3 x 10 8 A4q. This mass limit is 
comparable to the lowest mass galaxies in the ACSVCS. 
From this simulation, we selected a sample of low mass 
early- type cluster galaxies and their progenitors. We 
identified 126 massive galaxy clusters using the same cri- 
teria as in De Lucia et al. (2007), requiring a halo mass 
greater than 7 x 1§ 1a Mq at z = 0. This mass limit 
also corresponds roughly to the mass of Virgo (Bohringcr 
et al. 1994). Within each of these clusters, we then se- 
lected all galaxies with M z , s d ss > —19 that have early- 
type morphologies, matching the low luminosity sample 
we have in the ACSVCS. Early-type galaxies were defined 
by the bulge- to- total luminosity ratio with AM < 1.56 
(AM = Mb u ig e — Mtotai), using the empirical criteria of 
Simicn & de Vaucouleurs (1986). In addition, we selected 
only galaxies with g s dss — z sc i ss > 0.5 at z = in order to 
best match the colors of the ACSVCS galaxies. In total, 
we selected 15,506 simulated galaxies at z = and all of 
their progenitors over 63 snapshots back to z ~ 12. 

In Figure 18, we show the mass- weighted age of the 
stellar populations in the selected early-type dwarfs as a 
function of distance from the center of the cluster's dark 
matter halo. The error bars show the la width of the 
age distribution in each bin of radius. This figure clearly 
shows that the stellar age of the simulated dwarfs have 
a strong dependence on clustercentric radius, at least 
within a distance of 1.5 Mpc, with the central dwarfs 
having mean ages of 11.5 Gyr decreasing to 9 Gyr at the 
cluster outskirts. The radius within which the age gra- 
dient is evident is also the radius in Virgo within which 
we see elevated GC fractions. 

The central dwarfs in the simulation may be older, but 
they also need to form their stars more intensely (higher 
peak star formation rates and densities) in order to have 
more of their stars in massive clusters. Figure 19a shows 
the average normalized star formation rate of early-type 
cluster dwarfs as a function of lookback time. The star 
formation rate is normalized to the final stellar mass of 
the galaxy. The two curves represent SFR histories for 
dwarfs divided into two bins of clustercentric radius at 
1 Mpc from the cluster center. We note that these are 
average star formation histories, and that for any individ- 
ual galaxy, the bursts of star formation are more intense, 
brief, and stochastic. Combining galaxies allows us to see 
that the central dwarfs have a more peaked SFR at ear- 
lier times with a rapid falloff, whereas the outer dwarfs 
not only have a lower peak but also more star formation 
extending to later times. 

It has been suggested that massive star clusters form 
preferentially in high-pressure environments (Harris & 
Pudritz 1994; Elmegreen & Efremov 1997; Ashman & 
Zepf 2001), of which SFR surface density is one possi- 
ble indication. Figure 19b shows the ratio of SFR sur- 
face density in central dwarfs to that in outer dwarfs. 
We calculate the SFR surface density using the SFR and 
the disk radius, which the semi-analytic models calculate 
using the analytic model of Mo, Mao & White (1998). 
Not only is the peak star formation rate higher in cen- 
tral dwarfs, but the intensity of star formation, as mea- 
sured by the SFR surface density, is also higher during 
the epoch when these galaxies are forming most of their 
stars. At later times, the SFR surface density in the outer 
dwarfs is higher on average than that in central dwarfs, 



18 



Peng et al. 




4 6 8 10 
Lookback time (Gyr) 



14 



Fig. 19. — (a, top) Average normalized star formation rate versus 
lookback time for low mass early-type cluster galaxies in the Mil- 
lennium simulation. The solid line shows the average SFR history 
for "central" galaxies within 1 Mpc of the cluster center, and the 
dashed line shows "outer" galaxies. Low mass early-type galax- 
ies in the inner cluster regions formed their stars earlier and with 
higher peak star formation rate (SFR) than those in the outskirts. 
Higher peak SFRs could result in higher globular cluster mass frac- 
tions, (b, bottom) Average ratio in SFR surface density between 
inner and outer low mass early-type cluster galaxies as a func- 
tion of lookback time. Galaxies in the inner regions have a higher 
SFR surface density during the period of highest total star forma- 
tion (tiookback > 8 Gyr). This results in central galaxies having a 
higher mass fraction in massive star clusters. 

but this is at much lower absolute star formation rates 
and SFR surface densities. If GCs are formed in the same 
events that produced the bulk of these low mass galax- 
ies, then we would expect that the central dwarfs would 
retain a higher fraction of their stellar mass in massive 
star clusters. 

5.2.2. Inferring Cluster Formation Histories 

We can quantify these effects more directly by infer- 
ring the cluster formation rate (CFR) history, using an 
empirical relationship between Sm , the fraction of stellar 
mass formed in massive star clusters, and the SFR sur- 
face density, ~Esfr- Fitting the data of Larsen & Richtlcr 
(2000), we adopt Sm oc (^sfr) ' 8 , up to a maximum ef- 
ficiency of 100%. The cluster formation rate then scales 
with a combination of the SFR and the SFR surface den- 
sity as CFR oc SFR x (S S Ffl) ' 8 - The Larsen & Richtler 
(2000) data measured the fraction of luminosity in young 
massive clusters (YMCs), which had ages between 10 
and 500 Myr and more massive than ~ 3 x 10 A4®. 
Given that most clusters will be destroyed over a Hubble 
time, we need to estimate what fraction of these clus- 
ters will survive. We adopt the two-stage disruption 
law of Whitmore, Chandar & Fall (2007) where "infant 
mortality" causes clusters to disrupt in constant num- 
bers, dN/dr oc r -1 , for ages less than a few xlO 8 Gyr 
(Fall, Chandar & Whitmore 2005), and a constant mass 
loss over a Hubble time of fj, ev = 1.9 x 10 4 M ( 7 ) Gyr~ 1 
due to two-body relaxation (Fall & Zhang 2001; Jordan 
et al. 2007). Starting with the data of Larsen & Richtlcr 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



o 



Star Formation 
Cluster Formation 
SFR surface density 



R < 1 Mpc 




10 2 
10° 
10~ 2 

IO" 4 

10~ 6 
10° 
10' 

io- 2 
io- 3 

10~ 4 

IO 5 
W 6 



4 6 
Redshift 



10 



Fig. 20. — For central (top) and outer (bottom) early-type dwarfs 
in the Millennium Simulation, a comparison of the average star for- 
mation rate (solid), the average star formation rate surface density 
(dotted), and the inferred average cluster formation rate (dashed). 
The highest SFR surface densities occur at earlier times than the 
highest SFRs because the disks within which stars form arc smaller 
at high redshift. If massive star clusters preferentially form at high 
SFR surface densities, then GCs in early-type dwarfs will, on av- 
erage, have older ages and lower metallicities than the field stars. 
Formation rates are normalized to their maximum values. 

(2000), we estimate that ~ 1.8% of their young massive 
clusters survive to the present day. Using this formalism, 
we produce the corresponding cluster formation rate for 
central early- type dwarfs, normalized to the final mass 
in GCs. 

Figure 20 compares the average SFR and SFR surface 
density as a function of redshift for the sample of cen- 
tral dwarfs (top) and outer dwarfs (bottom). While the 
SFR peaks at z = 3.5 — 4, Y^sfr is highest at z = 10 
and falls a factor of 35 by z = 4. Because the formation 
rate of massive star clusters depends on both quantities, 
their peak epoch of formation will be earlier than that 
of the stars. The result is that the CFR peaks at earlier 
redshift, z ~ 4.5 to 5, than the total star formation rate, 
a difference of ~ 350-500 Myr. In all dwarfs, the de- 
pendence of massive star cluster formation on the SFR 
surface density naturally produces GCs that are older 
and more metal-poor than the stars. 

Depending on how massive star cluster formation 
scales with SFR and SFR density, we also expect that 
the mean age difference between the GCs and the field 
would be different in the central and outer dwarfs. Be- 
cause GCs in these dwarfs will mostly form at or before 
the SFR peak, they will always be old, but the field stars 
can continue to form at later times. This is evident when 
comparing the difference in lookback time for the SFR 
peaks for the central and outer simulated dwarfs in Fig- 
ure 19a and the difference in their mean ages in Figure 18. 
The outer dwarfs have mass- weighted mean ages that are 
younger by ~ 2.5 Gyr, but the times of their peak star 
formation rates differ by only 0.1 Gyr. Figure 20 shows 
how the CFR peak happens at similar redshifts for all 
dwarfs, but the peak in star formation happens later in 



Formation Efficiencies of Globular Clusters 



19 



CO 




0.0 0.5 1.0 1.5 2.0 2.5 

r (Mpc) 



Fig. 21. — Simulated Sm versus clustercentric radius for early- 
type dwarfs in the Millennium simulation. Colors from purple to 
red (or black to white) symbolize lower to higher numbers of simu- 
lated dwarfs in a particular cell in this diagram. Like in the obser- 
vations (Figures 10 and 11), there is a tendency for dwarf galaxies 
with higher Sm to be at the center of the cluster. 

outer dwarfs. One prediction is then that, on average, 
the difference between the mean age of the metal-poor 
GCs and the mean age of the stars should be larger in 
dwarf galaxies with lower GC mass fractions (i.e., the 
mean age of the galaxy is proportional to Sm)- 

Lastly, we test whether this framework can produce 
higher GC mass fraction in the central dwarfs. Fig- 
ure 21 shows Sm against clustercentric radius for sim- 
ulated dwarfs. The absolute scale for Sm is highly sensi- 
tive to the assumed destruction but is not important for 
our purposes. It is the relative comparison between cen- 
tral and outer simulated dwarfs that is of interest, and 
the high-SV; simulated dwarfs do appear preferentially 
in the central regions, similar to what the data shows in 
Figures 10 and 11. 

The simulations do not match the data exactly, as 
there also appears to be many Iow-Sm dEs within 1 Mpc, 
and the metallicities of the GCs as derived from the sim- 
ulation are generally too high, but it is encouraging that 
this simple scaling for cluster formation in the Millen- 
nium Simulation can reproduce many of the observed 
trends. A more detailed treatment of ram pressure and 
tidal stripping in cluster cores in a higher resolution sim- 
ulation may help resolve some discrepancies. 

5.3. Possible Mechanisms for Quenching Star 
Formation 

The central dwarfs may produce a higher fraction of 
GCs at or before their peak SFR, but it is also important 
that their subsequent star formation is rapidly quenched 
or kept at a low level. At these and lower masses, pho- 
toionization heating of the ISM by the UV background 
is believed to be an important mechanism for suppress- 
ing star formation in dwarf galaxies (Bullock, Kravtsov 
& Weinberg 2000; Benson et al. 2002), and this could 
universally truncate the epoch of efficient formation of 



metal-poor GCs, although detailed work still needs to be 
done to show that reionization is a plausible mechanism 
for halting GC formation. At least one scenario (Cen 
2001) has reionization as the cause for GC formation, an- 
other suggests that starburst-driven shocks may trigger 
GC formation (Scannapieco, Weisheit & Harlow 2004), 
and yet another shows that GCs themselves could plau- 
sibly be the source of reionizing photons (Ricotti 2004). 
So, the relationship between GC formation and reioniza- 
tion, if any, is still very much undetermined. 

One potential problem with the scenario where reion- 
ization halts GC formation is that in the simulations of 
Bekki et al. (2006), the bias introduced by the trunca- 
tion redshift is not very strong. Changing the truncation 
redshift from z = 15 to z = 6, for example, does not alter 
the Sn of galaxies by very much, and at least in these 
simulations, rcionization-induced biased GC formation 
by itself does not seem to be able to explain the scatter 
in Sn for dEs. However, it may be sufficient to explain 
the general trend towards higher GC fractions at lower 
masses seen in Figure 17. 

Another plausible explanation is that the environment 
plays an important role in quenching the later, lower level 
star formation that builds up the field. Ram pressure 
stripping of gas will be more efficient in low mass galaxies 
at the center of the cluster. Since these halos are the ear- 
liest to fall into the cluster, they will have their star for- 
mation quenched earlier and more efficiently than their 
counterparts in the outer cluster regions. Moore, Lake 
and Katz (1998) also showed through numerical simu- 
lations that "galaxy harassment" in clusters is efficient 
at transforming small spirals into low mass spheroidal 
galaxies. More recently, Mayer et al. (2007) used hydro- 
dynamic simulations to demonstrate that the environ- 
mental effects of tidal shocks and ram pressure stripping 
can combine to create the most dark matter-dominated 
dwarf spheroidals (dSphs) in the Local Group. Although 
the Virgo dEs are much more massive than dSphs such 
as Draco or Ursa Minor, it is possible that similar pro- 
cesses caused them to fail in their conversion of gas to 
field stars. Ram pressure stripping in Virgo has been ob- 
served out to a cluster radius of 0.8 Mpc (Kenney, van 
Gorkom & Vollmer 2004), similar to the radius within 
which we see enhanced GC fractions. The comparison of 
central dEs and those on the outskirts, however, shows 
no strong differences between their observed properties. 
If harassment was important, we might expect different 
morphologies or surface brightness profiles (Sersic n or 
r e ) in the two groups, but this is not the case. The cu- 
mulative GC color distributions of these two populations 
of dEs are also indistinguishable. If the central dEs were 
once much more luminous, we would expect them to be 
outliers on the color-magnitude relation, but the ones 
with high Sn are not. Only the innermost galaxies that 
are well within the M87 GC system and have no GCs are 
outliers and are red for their luminosities. These galaxies 
are good candidates for having been harassed. 

5.4. Globular Cluster Destruction 

Up until now, we have only discussed the possible 
variation in formation histories to explain the observed 
trends in GC fraction. However, star clusters can also be 
destroyed, and we expect that a substantial fraction of 
clusters initially formed will not survive a Hubble time. 
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Observations of young star clusters in nearby galaxies 
show that most star clusters are disrupted very early 
(Fall et al. 2005; Bastian et al. 2005). In addition, sub- 
sequent evolution through two-body relaxation, stellar 
mass loss, and tidal shocks will destroy even more low 
mass clusters (Fall & Rees 1977; Vesperini 1998; Fall & 
Zhang 2001). Globular clusters are thus the survivors 
of what was once a much larger population of star clus- 
ters. Any differences in the survival rates between differ- 
ent galaxies could produce different GC fractions in the 
present day. 

Unfortunately, the cosmological simulations used to 
study galaxy formation cannot at the present time model 
GC destruction. One benefit of studying trends in terms 
of Smi ° r GC mass fraction, is that while most of the 
destroyed star clusters after the initial 1 Gyr are pref- 
erentially low mass, most of the mass remains in the 
higher mass objects. Thus, GC mass fraction is more 
robust against disruption processes that preferentially 
destroy low mass clusters. If, however, there was a mech- 
anism that could affect the survival efficiency of massive 
clusters — perhaps variable infant mortality or dynamical 
friction as influenced by the dark matter halo profile — 
then destruction could play a role in driving the trends 
in GC mass fraction. This mechanism would require cen- 
tral dwarfs to be able to retain a larger fraction of the 
massive star clusters than their counterparts at larger 
cluster radii. This is perhaps counter to what one might 
expect if the cluster tidal field had a role in stripping or 
destroying GCs. We also see no signs in our data that 
central dwarfs have higher mean GC masses or different 
GC luminosity functions, things that might point to in- 
ternal destruction being a dominant mechanism driving 
GC fractions. 

For a real test of the mechanisms that drive GC forma- 
tion efficiency, it would be extremely beneficial to have 
more and better simulations of galaxies and their star 
cluster populations. Both hydrodynamic simulations and 
dark matter simulations coupled with semi-analytic mod- 
els can provide useful quantitative predictions. These 
would ideally be able to produce az = volume com- 
parable to the Virgo cluster with a stellar mass resolu- 
tion that would allow the resolution of GC mass objects 
(~ 10 5 X o ). 

5.5. Sm Variation in Metal-poor and Metal-rich GCs 

Another result of this study is the behavior of the spe- 
cific frequencies of blue and red GCs. The variation in 
Sm is dominated by the blue, metal-poor GCs, although 
the red GCs do exhibit similar trends with fewer GCs 
and lower significance. Figure 8 shows that the fraction 
of red GCs rarely exceeds 50% and that the trend of in- 
creasing red fraction with galaxy mass either flattens or 
turns over for the most massive galaxies. This suggests 
that massive ellipticals have accreted the most dwarf- 
like galaxies, the presumed original hosts of metal-poor 
GCs, or that the dwarfs which they accreted had the 
highest Sn, or both. Also, the rise in red GC fraction 
from M z = —19 to M z = —22 requires that SN, z ,biue and 
SN,z,red respectively fall and rise in equal proportions to 
produce the nearly constant Sn, z over this range. The 
fraction of blue GCs is lowest where the total Sn, z is 
at a minimum, implying again that they are driving the 
overall trends. 



The higher values of Spf,z,red for the massive galaxies, 
in particular M87, means that the formation of metal- 
rich GCs also occurs at different efficiencies. For massive 
galaxies, these results are consistent with the data of 
Rhode et al. (2005), who found that the mass- normalized 
number of GCs, T re d, is positively correlated with galaxy 
mass (shown in Figure 4 of Brodie & Strader 2006). 
Could these metal-rich GCs have been accreted from 
dwarfs in a scenario similar to that described for the 
metal-poor GCs? Figure 15 shows that the number frac- 
tion of red GCs may also be enhanced in central dwarfs, 
like the metal-poor GCs, but it is unclear whether a bi- 
asing scenario similar to the one described above provide 
an explanation since we would expect that the forma- 
tion of metal-rich GCs at a later time might actually 
be suppressed. However, the numbers we observe are 
very small and uncertain and may be sensitive to the 
dividing color between blue and red GCs, which is cho- 
sen to be the same in all galaxies at these luminosities. 
Another possibility is that the trend in the red GCs is 
determined more by the mass of the host, with red GCs 
forming in the most massive progenitor, and the more 
massive protogalaxics forming their stars earlier and with 
higher star formation rates. This is consistent with ob- 
servations of the stellar populations of elliptical galaxies, 
which find that the stars in massive ellipticals (L > L*) 
form have old ages and high [a/Fe] — i.e., they formed 
early and rapidly — while lower mass early-type galaxies 
have younger mean ages and lower [a/Fe] (Thomas et al. 
2005). Similar trends in [a/Fe] are seen in the GCs of 
these galaxies (Puzia et al. 2006). 

The fact that we see the same trends for both the blue 
and red GC sub-populations suggests the possibility that 
the two populations arc not very distinct, and that their 
separation is merely the result of placing a dividing line 
on a continuum of GC properties. This possibility has 
been raised by Yoon, Yi, & Lee (2006), who suggest that 
the GC mctallicity bimodality in early-type galaxies is a 
result of a nonlinear mctallicity-color relation. Although 
our data does not constrain this hypothesis, we do not 
observe a clear distinction between blue and red GC frac- 
tions across galaxy mass or environment. 

6. CONCLUSIONS 

We have measured globular cluster specific frequencies 
(Sn), luminosity fractions (fix), and stellar mass frac- 
tions (Sm) f° r 100 early- type galaxies in the ACS Virgo 
Cluster Survey. These galaxies span the mass range from 
giants to dwarfs (—22 < Mb < —15) and these represent 
the largest homogeneous catalog of GC number and mass 
fractions to date. We have studied these quantities as a 
function of galaxy mass and environment and find that: 

1. Globular cluster fractions can be high (Sjv, z > 2) 
for both high and low luminosity early-type galax- 
ies, but are universally low for intermediate lumi- 
nosity galaxies (-22 < M z < -19). 

2. There is a large spread in GC fraction in early- type 
dwarfs (0 < Sn.z < 5.5) which can be understood 
as an underlying dependence on environment. Al- 
most all dwarfs (M z > —19) with Sn, z > 2 are 
within a projected radius of 1 Mpc from M87 and 
the cluster center (or 1.5 Mpc in three-dimensional 
cluster radius). The spatial distributions of GCs 
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in the high-Sjv dwarfs are centrally concentrated, 
showing that they are intrinsic to the host dwarf 
and not interloping GCs from M87 or an intra- 
cluster population. We do not detect higher GC 
fractions in dwarfs around the other massive ellip- 
ticals in the cluster. We present this as evidence 
that GC formation in low mass galaxies is biased 
towards the densest environments. 

3. Galaxies within R p ~ 40 kpc of M87 and M49 have 
few or no GCs, and are likely to have had their 
GC systems tidally stripped by their giant neigh- 
bors. Galaxies out to R p ~ 100 kpc from M87 
(VCC 1185) also appear to be affected. 

4. Analyzing the blue and red GC populations sep- 
arately, we find that the fraction of red GCs in- 
creases with galaxy luminosity until M z ~ —22 at 
which point it flattens or declines. Trends in GC 
fraction are dominated by the blue GCs, although 
the iSjv,z of red GCs does exhibit similar relative en- 
hancement for massive galaxies, and possibly shows 
a weak enhancement for central dwarfs. 

5. We use a globally averaged Mh/ Mg*~Mg* rela- 
tion from halo occupation studies to test whether 
trends in GC fraction can be explained by the as- 
sumption that the mass in GCs is directly propor- 
tional to the halo mass, Aiac -M-h- While this 
may be able to explain the mean trend in Sm-i it 
is unable to account simultaneously for all dwarfs, 
particularly the dwarfs with high GC fractions in 
the cluster's central regions. 

6. Comparisons with semi-analytic models of galaxy 
formation in the Millennium Simulation show that 
early-type dwarfs in the central 1 Mpc of massive 
galaxy clusters are expected to be older and have 
higher peak star formation rates and SFR surface 
densities than their counterparts on the outskirts 
of the cluster. The higher stellar mass fractions in 
globular cluster for central dwarfs can be explained 
if higher SFR surface densities are responsible for 
more efficient formation of massive star clusters. 

7. The peak SFR surface density in simulated dwarfs 
occurs before the peak SFR, which we propose as 
an explanation for why GCs are, on average, both 
older and more metal-poor than the field stars in 
their host galaxies. 

We present a picture of globular cluster system for- 
mation in the Virgo Cluster where the highest GC mass 
fractions are formed in the oldest systems. These pro- 
genitors formed a larger fraction of their stars at higher 
peak SFR surface densities, and also have star formation 
suppressed at later times. In regions of high density, 
halos collapse and star formation starts earlier, and the 
mechanisms for truncating and suppressing subsequent 
star formation are stronger. For all dwarfs, the highest 
SFR surface densities occur earlier than the peak SFR, 
naturally producing the metal-poor GC populations we 
see today — old (> 10 Gyr) and more metal-poor than the 
bulk of the field stars. Merging of many low mass progen- 
itors with high GC fractions produce the extremely high 
Sn for blue GCs seen in the cD galaxy, M87. They may 



also, at a lower level, produce the elevated blue GC frac- 
tions seen in other massive ellipticals. Future detailed 
simulations of GC system formation will be crucial to 
test, in a quantitative way, the scenarios of galaxy and 
GC system formation that we are beginning to assemble. 

E. W. P. gratefully acknowledges the support of the 
National Research Council of Canada's Plaskett Re- 
search Fellowship at the Herzberg Institute of Astro- 
physics. We thank Gabriella De Lucia, Gerard Lemson, 
and Darren Croton for their help with the Millennium 
Simulation, and thank Rupali Chandar, Dean McLaugh- 
lin, Anil Seth, and Peter Stetson for useful discussions 
and comments on the manuscript. We also thank Pe- 
ter Stetson for providing point spread functions used in 
the analysis of our WFPC2 data. We thank Naoyuki 
Tamura, who provided the numbers for their M87 GC 
radial profile. We are grateful to the anonymous referee 
who provided comments that improved the manuscript. 
M. J. W. and M. T. acknowledge support through NSF 
grant AST 02-05960. 

This publication makes use of data products from the 
Sloan Digital Sky Survey (SDSS). Funding for SDSS 
and SDSS-II has been provided by the Alfred P. Sloan 
Foundation, the Participating Institutions, the National 
Science Foundation, the U.S. Department of Energy, 
the National Aeronautics and Space Administration, the 
Japanese Monbukagakusho, the Max Planck Society, and 
the Higher Education Funding Council for England. The 
SDSS Web site is http://www.sdss.org/. 

The SDSS is managed by the Astrophysical Research 
Consortium (ARC) for the Participating Institutions. 
The Participating Institutions are the American Mu- 
seum of Natural History, Astrophysical Institute Pots- 
dam, University of Basel, University of Cambridge, Case 
Western Reserve University, The University of Chicago, 
Drexel University, Fermilab, the Institute for Advanced 
Study, the Japan Participation Group, The Johns Hop- 
kins University, the Joint Institute for Nuclear Astro- 
physics, the Kavli Institute for Particle Astrophysics 
and Cosmology, the Korean Scientist Group, the Chi- 
nese Academy of Sciences (LAMOST), Los Alamos Na- 
tional Laboratory, the Max-Planck-Institute for Astron- 
omy (MPIA) , the Max-Planck-Institute for Astrophysics 
(MPA), New Mexico State University, Ohio State Univer- 
sity. University of Pittsburgh, University of Portsmouth, 
Princeton University, the United States Naval Observa- 
tory, and the University of Washington. 

This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 

This research has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronautics and 
Space Administration. 

Facilities: HST(ACS,WFPC2) 



22 



Peng et al. 



REFERENCES 



Adelman-McCarthy, J. K., et al. 2007, ApJS, 172, 634 

Ashman, K. M., & Zepf, S. E. 1992, ApJ, 384, 50 

Ashman, K. M., & Zepf, S. E. 1998, Globular cluster systems 

/ Keith M. Ashman, Stephen E. Zepf. Cambridge, U. K. ; 

New York : Cambridge University Press, 1998. (Cambridge 

astrophysics series ; 30) QB853.5 .A84 1998 
Ashman, K. M., & Zepf, S. E. 2001, AJ, 122, 1888 
Baggct, S. et al. 2002, in HST WFPC2 Data Handbook, v. 4.0, ed. 

B. Mobasher, Baltimore, STScI 
Bastian, N., Gieles, M., Lamcrs, H. J. G. L. M., Scheepmakcr, 

R. A., & de Grijs, R. 2005, A&A, 431, 905 
Beasley, M. A., Baugh, C. M., Forbes, D. A., Sharpies, R. M., & 

Frenk, C. S. 2002, MNRAS, 333, 383 
Bckki, K., Yahagi, H., & Forbes, D. A. 2006, ApJ, 645, L29 
Benson, A. J., Frenk, C. S., Lacey, C. G., Baugh, C. M., & Cole, 

S. 2002, MNRAS, 333, 177 
Berlind, A. A., & Weinberg, D. H. 2002, ApJ, 575, 587 
Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Binggeli, B., Sandage, A., & Tammann, G. A. 1985, AJ, 90, 1681 
Blakcslee, J. P. 1997, ApJ, 481, L59 

Blakeslee, J. P., Tonry, J. L., & Metzger, M. R. 1997, AJ, 114, 482 
Blakeslee, J. P. 1999, AJ, 118, 1506 

Blakeslee, J. P., Meurer, G. R., Ford, H. C, Benitez, N., White, 
R. L., Zekser, K. C, & Sirianni, M. 2003, Extragalactic Globular 
Clusters and their Host Galaxies, 25th meeting of the IAU, Joint 
Discussion 6, 17 July 2003, Sydney, Australia, 6 

Blakeslee, J. P., Anderson, K. R., Meurer, G. R., Bcm'tez, N., 
& Magee, D. 2003, Astronomical Data Analysis Software and 
Systems XII, 295, 257 

Blanton, M. R., & Roweis, S. 2007, AJ, 133, 734 

Bohringer, H., Briel, U. G., Schwarz, R. A., Voges, W., Hartner, 
G., & Trumper, J. 1994, Nature, 368, 828 

Brodie, J. P., & Huchra, J. P. 1991, ApJ, 379, 157 

Brodie, J. P., & Strader, J. 2006, ARA&A, 44, 193 

Bullock, J. S., Kravtsov, A. V., & Weinberg, D. H. 2000, ApJ, 539, 
517 

Cantiello, M., & Blakcslee, J. P. 2007, ApJ, 669, 982 

Cen, R. 2001, ApJ, 560, 592 

Cappellari, M., et al. 2006, MNRAS, 366, 1126 

Chaboyer, B., Demarque, P., Kernan, P. J., & Krauss, L. M. 1996, 

Science, 271, 957 
Chabrier, G. 2003, PASP, 115, 763 

Chandar, R., Whitmore, B., & Lee, M. G. 2004, ApJ, 611, 220 
Cote, P., Marzke, R. O., & West, M. J. 1998, ApJ, 501, 554 
Cote, P., Marzke, R. O., West, M. J., & Minniti, D. 2000, ApJ, 
533, 869 

Cote, P., et al. 2001, ApJ, 559, 828 

Cote, P., McLaughlin, D. E., Cohen, J. G., & Blakeslee, J. P. 2003, 

ApJ, 591, 850 
Cote, P., et al. 2004, ApJS, 153, 223 (Paper I) 
Cote, P., ct al. 2006, ApJS, 165, 57 
Croton, D. J., ct al. 2006, MNRAS, 365, 11 
Cote, P., et al. 2007, ApJ, submitted 

De Lucia, G., Springel, V., White, S. D. M., Croton, D., & 

Kauffmann, G. 2006, MNRAS, 366, 499 
De Lucia, G., & Blaizot, J. 2007, MNRAS, 375, 2 
de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Jr., Buta, 

R. J., Paturel, G., & Fouque, P. 1991, Volume 1-3, XII, 2069 

pp. 7 figs.. Springer- Verlag Berlin Heidelberg New York 
Dekel, A., & Birnboim, Y. 2006, MNRAS, 368, 2 
Diemand, J., Madau, P., & Moore, B. 2005, MNRAS, 364, 367 
Dirsch, B., Richtler, T., Geisler, D., Forte, J. C, Bassino, L. P., & 

Gieren, W. P. 2003a, AJ, 125, 1908 
Dirsch, B., Richtler, T., & Bassino, L. P. 2003b, A&A, 408, 929 
Dirsch, B., Schuberth, Y., & Richtler, T. 2005, A&A, 433, 43 
Durrcll, P. R., Harris, W. E., Geisler, D., & Pudritz, R. E. 1996a, 

AJ, 112, 972 

Durrell, P. R., McLaughlin, D. E., Harris, W. E., & Hanes, D. A. 

1996b, ApJ, 463, 543 
Elmegreen, B. G., & Efremov, Y. N. 1997, ApJ, 480, 235 
Fall, S. M., & Rees, M. J. 1977, MNRAS, 181, 37P 
Fall, S. M., & Zhang, Q. 2001, ApJ, 561, 751 

Fall, S. M., Chandar, R., & Whitmore, B. C. 2005, ApJ, 631, L133 
Ferrarese, L., et al. 2006a, ApJS, 164, 334 



Ferrarese, L., et al. 2006b, ApJ, 644, L21 

Forbes, D. A., Brodie, J. P., & Huchra, J. 1996, AJ, 112, 2448 

Forbes, D. A. 2005, ApJ, 635, L137 

Ford, H. C, et al. 1998, Proc. SPIE, 3356, 234 

Freedman, W. L., et al. 2001, ApJ, 553, 47 

Gao, L., Springel, V., & White, S. D. M. 2005, MNRAS, 363, L66 
Geha, M., Guhathakurta, P., & van der Marel, R. P. 2003, AJ, 126, 
1794 

Goudfrooij, P., Strader, J., Brcnncman, L., Kissler-Patig, M., 

Minniti, D., & Edwin Huizinga, J. 2003, MNRAS, 343, 665 
Gomez, M., & Richtler, T. 2004, A&A, 415, 499 
Hasxgan, M., et al. 2005, ApJ, 627, 203 (Paper VII) 
Harris, G. L. H., Harris, W. E., & Geisler, D. 2004, AJ, 128, 723 
Harris, W. E., & van den Bergh, S. 1981, AJ, 86, 1627 
Harris, W. E. 1991, ARA&A, 29, 543 
Harris, W. E., & Pudritz, R. E. 1994, ApJ, 429, 177 
Harris, W. E., Harris, G. L. H., & McLaughlin, D. E. 1998, AJ, 
115, 1801 

Harris, W. E., & Harris, G. L. H. 2002, AJ, 123, 3108 

Hoekstra, H., Hsieh, B. C, Yee, H. K. C, Lin, H., & Gladders, 

M. D. 2005, ApJ, 635, 73 
Jarrett, T. H., Chester, T., Cutri, R., Schneider, S. E., & Huchra, 

J. P. 2003, AJ, 125, 525 
Jordan, A., West, M. J., Cote, P., & Marzke, R. O. 2003, AJ, 125, 

1642 

Jordan, A., Cote, P., West, M. J., Marzke, R. O., Minniti, D., & 

Rejkuba, M. 2004a, AJ, 127, 24 
Jordan, A., et al. 2004a, ApJS, 154, 509 (Paper II) 
Jordan, A., ct al. 2004b, ApJ, 613, 279 (Paper III) 
Jordan, A., ct al. 2005, ApJ, 634, 1002 (Paper X) 
Jordan, A., et al. 2006, ApJ, 651, L25 
Jordan, A., ct al. 2007, ApJS, 171, 101 (Paper XII) 
Kauffmann, G., et al. 2003, MNRAS, 341, 54 

Kenney, J. D. P., van Gorkom, J. H., & Vollmer, B. 2004, AJ, 127, 
3361 

King, I. R. 1966, AJ, 71, 64 

Kodama, T., Arimoto, N., Barger, A. J., & Arag'on-Salamanca, A. 

1998, A&A, 334, 99 
Koekemoer, A. M., Fruchtcr, A. S., Hook, R. N, & Hack, 
W. 2002, The 2002 HST Calibration Workshop : Hubble 
after the Installation of the ACS and the NICMOS Cooling 
System, Proceedings of a Workshop held at the Space Telescope 
Science Institute, Baltimore, Maryland, October 17 and 18, 
2002. Edited by Santiago Arribas, Anton Koekemoer, and Brad 
Whitmore. Baltimore, MD: Space Telescope Science Institute, 
2002., p.339, 339 
Kundu, A., & Whitmore, B. C. 2001a, AJ, 121, 2950 
Kundu, A., & Whitmore, B. C. 2001b, AJ, 122, 1251 
Kravtsov, A. V., & Gncdin, O. Y. 2005, ApJ, 623, 650 
Kronawitter, A., Saglia, R. P., Gerhard, O., & Bender, R. 2000, 

A&AS, 144, 53 
Larsen, S. S., & Richtler, T. 2000, A&A, 354, 836 
Larsen, S. S., Brodie, J. P., Huchra, J. P., Forbes, D. A., & 

Grillmair, C. J. 2001, AJ, 121, 2974 
Lotz, J. M., Miller, B. W., & Ferguson, H. C. 2004, ApJ, 613, 262 
Mandclbaum, R., Seljak, U., Kauffmann, G., Hirata, C. M., & 

Brinkmann, J. 2006, MNRAS, 368, 715 
Mayer, L., Kazantzidis, S., Mastropictro, C, & Wadsley, J. 2007, 

Nature, 445, 738 
McLaughlin, D. E. 1999a, AJ, 117, 2398 
McLaughlin, D. E. 1999b, ApJ, 512, L9 
Mei, S., et al. 2005a, ApJS, 156, 113 
Mei, S., et al. 2005b, ApJ, 625, 121 (Paper V) 
Mei, S., et al. 2007, ApJ, 655, 144 
Micske, S., et al. 2006, ApJ, 653, 193 (Paper XIV) 
Miller, B. W., Lotz, J. M., Ferguson, H. C, Stiavelli, M., & 

Whitmore, B. C. 1998, ApJ, 508, L133 
Miller, B. W., & Lotz, J. M. 2007, ApJ, 670, 1074 
Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 
Moore, B., Lake, G., & Katz, N. 1998, ApJ, 495, 139 
Moore, B., Diemand, J., Madau, P., Zemp, M., & Stadel, J. 2006, 

MNRAS, 368, 563 
Napolitano, N. R., et al. 2005, MNRAS, 357, 691 
Peng, E. W., Ford, H. C, & Freeman, K. C. 2004, ApJ, 602, 685 



Formation Efficiencies of Globular Clusters 



23 



Peng, E. W., et al. 2006a, ApJ, 639, 95 (Paper IX) 
Peng, E. W., et al. 2006b, ApJ, 639, 838 (Paper XI) 
Puzia, T. H., Kissler-Patig, M., & Goudfrooij, P. 2006, ApJ, 648, 
383 

Puzia, T. H., & Sharina, M. E. 2007, ArXiv e-prints, 710, 

arXiv:0710.1550 
Rhode, K. L., & Zepf, S. E. 2001, AJ, 121, 210 
Rhode, K. L., & Zepf, S. E. 2004, AJ, 127, 302 
Rhode, K. L., Zepf, S. E., & Santos, M. R. 2005, ApJ, 630, L21 
Ricotti, M. 2004, The Formation and Evolution of Massive Young 

Star Clusters, 322, 509 
Romanowsky, A. J., Douglas, N. G., Arnaboldi, M., Kuijken, K., 

Merrifield, M. R., Napolitano, N. R., Capaccioli, M., & Freeman, 

K. C. 2003, Science, 301, 1696 
Salpcter, E. E. 1955, ApJ, 121, 161 
Guzik, J., & Seljak, U. 2002, MNRAS, 335, 311 
Santos, M. R. 2003, Extragalactic Globular Cluster Systems, 348 
Schlegcl, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Searle, L., & Zinn, R. 1978, ApJ, 225, 357 

Sersic, J. L. 1968, Cordoba, Argentina: Observatorio Astronomico, 
1968 

Scth, A., Olscn, K., Miller, B., Lotz, J., & Telford, R. 2004, AJ, 
127, 798 

Sirianni, M., Jee, M.J., Bentez, N., Blakeslee, J. P., Martel, A.R., 
Meurer, G., Clampin, M., De Marchi, G., Ford, H.C., Gilliland, 
R., Hartig, G.F., Illingworth, G.D., Mack, J., & McCann, W.J. 
2005, PASP, accepted 

Sivakoff, G. R., et al. 2007, ApJ, 660, 1246 

Skrutskie, M. F., et al. 2006, AJ, 131, 1163 

Spergel, D. N., et al. 2007, ApJS, 170, 377 



Spitlcr, L. R., Forbes, D. A., Strader, J., Brodie, J. P., & Gallagher, 

J. S., Ill 2007, ArXiv e-prints, 712, arXiv:0712.1382 
Springel, V., et al. 2005, Nature, 435, 629 

Strader, J., Brodie, J. P., Spitler, L., & Beasley, M. A. 2006, AJ, 
132, 2333 

Thomas, D., Maraston, C, Bender, R., & Mcndes de Oliveira, C. 

2005, ApJ, 621, 673 
Tonry, J. L., Dressier, A., Blakeslee, J. P., Ajhar, E. A., Fletcher, 

A. B., Luppino, G. A., Metzger, M. R., & Moore, C. B. 2001, 

ApJ, 546, 681 
Tremonti, C. A., et al. 2004, ApJ, 613, 898 

Vale, A., & Ostriker, J. P. 2007, ArXiv Astrophysics e-prints, 

arXiv:astro-ph/0701096 
van den Bergh, S. 1975, ARA&A, 13, 217 

van den Bosch, F. C, Yang, X., & Mo, H. J. 2003, MNRAS, 340, 
771 

van Dokkum, P. G. 2001, PASP, 113, 1420 
Vesperini, E. 1998, MNRAS, 299, 1019 
West, M. J. 1993, MNRAS, 265, 755 

West, M. J., Cote, P., Jones, C, Forman, W., & Marzke, R. O. 

1995, ApJ, 453, L77 
West, M. J., Cote, P., Marzke, R. O., & Jordan, A. 2004, Nature, 

427, 31 

West, A. A., et al. 2007, AJ, submitted 

Whitmore, B. C, Chandar, R., & Fall, S. M. 2007, AJ, 133, 1067 
Williams, B. F., et al. 2007, ApJ, 654, 835 
Yoon, S.-J., Yi, S. K., & Lee, Y.-W. 2006, Science, 311, 1129 
Zepf, S. E., Ashman, K. M., & Geisler, D. 1995, ApJ, 443, 570 



Peng et al. 



TABLE l 

Global Properties of ACS Virgo Cluster Survey Galaxies 



No. 
(1) 


VCC 
(2) 


M v 
(3) 


M z 
(4) 


R 

(5) 


R'iD 

(6) 


(7) 


M* 
(8) 


1 


1226 


-22.90 


-23.87 


1.27 


1.28 ± 0.39 


263.04 


531.84 ± 110.48 


2 


1316 


-22.66 


-23.86 


0.00 


0.00 


188.40 


302.27 ±79.13 


3 


1978 


-22.41 


-23.42 


0.94 


1.04 ±0.46 


167.62 


339.32 ± 50.28 


4 


881 


-22.54 


-23.53 


0.36 


1.49 ± 0.42 


181.42 


289.63 ± 59.87 


5 


798 


-22.34 


-23.25 


1.71 


2.09 ± 0.42 


132.84 


186.50 ± 43.84 


6 


763 


-22.29 


-23.20 


0.43 


2.09 ± 0.44 


141.95 


236.49 ±61.04 


7 


731 


-22.31 


-23.37 


1.53 


6.76 ±0.54 


141.76 


226.21 ± 52.45 


8 


1535 


-21.38 


-22.40 


1.37 




57.88 


77.21 ± 16.21 


9 


1903 


-21.18 


-22.19 


0.82 


1.87 ±0.35 


50.26 


83.87 ± 19.10 


10 


1632 


-21.36 


-22.33 


0.34 


0.63 ± 0.38 


60.33 


95.35 ± 16.89 


11 


1231 


-20.70 


-21.58 


0.31 


1.30 ± 0.36 


32.37 


53.46 ± 12.30 


12 


2095 


-20.78 


-20.95 


1.59 




33.74 


52.75 ± 10.12 


13 


1154 


-20.85 


-21.79 


0.47 


0.66 ± 0.38 


36.86 


77.95 ± 13.64 


14 


1062 


-20.56 


-21.34 


0.77 


1.42 ± 0.36 


28.65 


52.69 ±8.88 


15 


2092 


-20.46 


-21.64 


1.54 


1.61 ± 0.38 


25.55 


47.32 ± 7.66 


16 


369 


-20.19 


-20.41 


0.79 


1.04 ±0.37 


20.17 


32.46 ± 7.67 


17 


759 


-20.44 


-21.41 


0.46 


0.60 ± 0.40 


25.74 


45.15 ± 10.04 


18 


1692 


-20.26 


-21.04 


1.54 


1.64 ±0.40 


21.85 


33.97 ±6.55 


19 


1030 


-20.45 


-21.38 


0.30 


0.33 ± 0.39 


22.63 


13.32 ±9.05 


20 


2000 


-19.60 


-20.66 


1.03 


2.03 ±0.35 


11.71 


23.76 ±4.10 


21 


685 


-20.11 


-21.04 


1.33 




19.12 


31.23 ±6.69 


22 


1664 


-19.90 


-20.97 


0.48 


0.77 ±0.38 


15.39 


26.00 ±6.15 


23 


654 


-19.88 


-20.59 


1.35 




14.22 


22.74 ±4.12 


24 


944 


-19.80 


-20.66 


0.86 


1.05 ± 0.38 


13.94 


29.14 ±4.88 


25 


1938 


-19.93 


-20.86 


0.89 


1.00 ±0.40 


15.71 


24.86 ±5.50 


26 


1279 


-19.78 


-20.72 


0.04 


0.47 ± 0.40 


13.42 


22.03 ± 4.03 


27 


1720 


-19.63 


-20.68 


0.94 


0.97 ± 0.38 


11.55 


20.21 ± 5.31 


28 


355 


-19.42 


-20.41 


1.07 


1.55 ±0.36 


10.06 


15.78 ±3.62 


29 


1619 


-19.36 


-20.19 


0.33 


1.51 ±0.36 


8.65 


17.24 ±3.12 


30 


1883 


-19.68 


-20.73 


1.65 


1.69 ±0.38 


11.19 


16.66 ± 3.92 


31 


1242 


-19.38 


-20.29 


0.49 


1.42 ±0.43 


9.07 


15.27 ±3.17 


32 


784 


-19.31 


-20.26 


1.00 


1.30 ±0.37 


8.48 


16.79 ±3.14 


33 


1537 


-18.99 


-19.95 


0.39 


1.18 ± 0.36 


6.18 


10.24 ± 2.35 


34 


778 


-19.56 


-20.03 


0.79 


1.17 ±0.41 


10.73 


18.18 ±4.19 


35 


1321 


-18.80 


-19.93 


1.26 


1.99 ±0.35 


4.86 


6.98 ± 1.75 


36 


828 


-19.15 


-20.07 


0.38 


1.42 ±0.42 


7.51 


13.69 ±3.00 


37 


1250 


-18.98 


-20.10 


0.06 


0.79 ± 0.57 


5.23 


3.68 ± 2.04 


38 


1630 


-19.07 


-20.06 


0.34 


0.51 ±0.38 


7.03 


11.61 ±2.39 


39 


1146 


-18.93 


-19.64 


0.28 


0.60 ± 0.38 


5.71 


8.69 ±2.17 


40 


1025 


-19.58 


-20.57 


1.24 


5.98 ±0.53 


10.90 


21.39 ±4.25 


41 


1303 


-18.84 


-19.80 


0.97 


0.98 ± 0.39 


5.55 


10.39 ± 2.33 


42 


1913 


-18.87 


-19.74 


1.58 


1.67 ±0.40 


5.67 


10.74 ±2.72 


43 


1327 


-19.10 


-19.81 


0.04 


1.45 ± 0.51 


8.36 


17.24 ±2.34 


44 


1125 


-19.07 


-19.44 


0.24 




6.80 


8.06 ± 3.40 


45 


1475 


-18.56 


-19.42 


1.13 


1.15 ± 0.38 


4.07 


7.73 ± 1.50 


46 


1178 


-18.35 


-19.35 


1.22 


1.48 ±0.37 


3.64 


7.09 ± 1.42 


47 


1283 


-18.65 


-19.70 


0.34 


0.86 ±0.41 


4.83 


9.08 ± 1.79 


48 


1261 


-18.42 


-19.21 


0.47 


1.44 ±0.51 


3.14 


4.87 ± 1.69 


49 


698 


-18.78 


-19.71 


0.58 


1.80 ±0.43 


5.03 


9.52 ±2.56 


50 


1422 


-17.97 


-18.73 


0.62 


1.60 ±0.43 


2.18 


3.82 ± 1.35 


51 


2048 


-17.85 


-18.57 


1.32 




1.89 


2.94 ± 1.08 


52 


1871 


-17.31 


-18.32 


0.79 


1.43 ± 0.44 


1.36 


2.26 ±0.58 


53 


9 


-18.04 


-18.80 


1.57 


1.64 ±0.65 


2.17 


3.08 ±0.61 


54 


575 


-18.42 


-19.33 


1.34 


5.30 ±0.61 


3.56 


5.00 ± 1.75 


55 


1910 


-17.39 


-18.33 


0.82 


1.11 ± 0.37 


1.43 


2.10 ±0.81 


56 


1049 


-16.69 


-17.34 


1.26 


1.43 ±0.53 


0.60 


0.52 ±0.40 


57 


856 


-17.57 


-18.25 


0.76 


0.76 ± 0.47 


1.47 


2.22 ±0.85 


58 


140 


-17.51 


-18.27 


1.24 


1.30 ±0.46 


1.37 


2.33 ±0.71 


59 


1355 


-17.51 


-18.16 


0.50 


0.52 ±0.63 


1.29 


1.82 ±0.35 


60 


1087 


-17.79 


-18.64 


0.25 


0.34 ± 0.46 


1.84 


3.29 ± 1.07 


61 


1297 


-17.67 


-18.75 


0.04 


0.23 ± 0.46 


1.95 


3.96 ±0.62 


62 


1861 


-17.60 


-18.34 


0.80 


1.15 ± 0.45 


1.63 


2.88 ± 1.00 


63 


543 


-17.41 


-18.19 


0.91 


1.44 ±0.44 


1.27 


2.19 ±0.65 


64 


1431 


-17.39 


-18.25 


0.34 


0.90 ± 0.45 


1.37 


2.20 ±0.74 


65 


1528 


-17.16 


-18.04 


0.34 


0.70 ± 0.45 


1.06 


1.63 ±0.48 


66 


1695 


-17.49 


-18.32 


0.43 


0.46 ± 0.54 


1.30 


1.69 ±0.78 


67 


1833 


-17.13 


-17.89 


1.22 


1.33 ±0.45 


0.97 


0.86 ±0.48 


68 


437 


-17.82 


-18.69 


1.62 


1.65 ± 0.48 


1.92 


2.80 ± 1.25 


69 


2019 


-17.36 


-18.18 


1.08 


1.13 ±0.48 


1.17 


1.02 ±0.72 


70 


33 


-16.39 


-17.01 


1.48 


2.19 ±0.57 


0.46 


0.43 ± 0.25 


71 


200 


-17.12 


-17.80 


1.02 


1.70 ±0.59 


0.95 


1.34 ±0.49 


72 


571 


-17.32 


-18.29 


1.40 


7.27 ± 1.01 


1.07 


0.98 ±0.65 


73 


21 


-16.83 


-17.38 


1.58 




0.58 


0.56 ±0.37 


74 


1488 


-16.78 


-17.34 


0.88 




0.57 


0.41 ± 0.38 


75 


1779 


-16.90 


-17.38 


0.89 




0.63 


1.02 ±0.38 
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TABLE 1 — Continued 



No. 
(1) 


VCC 
(2) 


M v 
(3) 


M z 
(4) 


R 

(5) 


R 3D 
(6) 


(7) 


M* 
(8) 


76 


1895 


-16.60 


-17.29 


1.16 


1.47 ±0.37 


0.57 


0.76 ± 0.39 


77 


1499 


-16.53 


-16.95 


0.22 




0.34 


0.14 ±0.19 


78 


1545 


-16.91 


-17.74 


0.26 


0.30 ± 0.54 


0.87 


1.41 ±0.42 


79 


1192 


-16.86 


-18.14 


1.27 




0.94 


1.85 ±0.62 


80 


1857 


-16.64 


-17.32 


0.90 




0.53 


0.73 ±0.14 


81 


1075 


-16.78 


-17.55 


0.63 


0.87 ± 0.68 


0.73 


1.11 ±0.47 


82 


1948 


-16.06 


-16.78 


0.99 




0.35 


0.50 ± 0.09 


83 


1627 


-16.46 


-17.39 


0.34 


2.19 ±0.34 


0.59 


1.03 ±0.32 


84 


1440 


-16.86 


-17.66 


0.88 


1.21 ±0.45 


0.78 


1.19 ±0.44 


85 


230 


-16.21 


-16.96 


0.96 


1.41 ±0.84 


0.43 


0.69 ±0.29 


86 


2050 


-16.36 


-17.15 


1.16 


1.47 ±0.52 


0.47 


0.29 ±0.26 


87 


1993 


-16.30 


-17.02 


0.95 


1.03 ±0.38 


0.48 


0.75 ±0.27 


88 


751 


-16.97 


-17.92 


1.72 


2.05 ± 0.44 


0.92 


1.41 ±0.51 


89 


1828 


-16.78 


-17.47 


0.67 


0.67 ±0.55 


0.72 


1.26 ±0.47 


90 


538 


-16.24 


-17.10 


1.62 


6.47 ±0.86 


0.44 


0.63 ±0.26 


91 


1407 


-16.72 


-17.43 


0.17 


0.23 ±0.39 


0.68 


1.24 ±0.42 


92 


1886 


-16.25 


-16.93 


0.76 




0.36 


0.51 ±0.10 


93 


1199 


-15.47 


-16.94 


1.25 




0.29 


0.58 ±0.16 


94 


1743 


-16.33 


-16.95 


0.84 


1.27 ±0.75 


0.45 


0.40 ± 0.29 


95 


1539 


-16.05 


-17.12 


0.25 


0.26 ±0.87 


0.37 


0.52 ±0.11 


96 


1185 


-16.77 


-17.37 


0.10 


0.10 ±0.79 


0.78 


0.89 ±0.51 


97 


1826 


-16.03 


-16.71 


0.98 


1.11 ±0.61 


0.35 


0.61 ±0.22 


98 


1512 


-16.25 


-16.54 


0.38 


1.50 ±0.42 


0.43 


0.61 ±0.12 


99 


1489 


-15.61 


-16.41 


0.45 




0.22 


0.32 ±0.06 


100 


1661 


-15.81 


-17.40 


0.70 


1.26 ± 1.20 


0.29 


0.43 ± 0.08 



1 Running number, sorted by increasing Bt magnitude 

2 Number in Virgo Cluster Catalog 

3 Absolute V magnitude 

4 Absolute z magnitude 

5 Projected distance from M87 (VCC 1316), in Mpc 

6 3-dimcnsional distance from M87 (VCC 1316), in Mpc, using polynomial calibration from Paper XIII 

7 Stellar z luminosity (10 9 Lq) 

8 Stellar mass (1O 9 M ) 



TABLE 2 

GC Formation Efficiencies in ACS Virco Cluster Survey Galaxies 



to 



No. 
(1) 



vcc 

(2) 



(3) 



(4) 



S JV,2 

(5) 



(6) 



(7) 



S M 
(S) 



O) 



S N,z,red 
(10) 



J 6(ue 
(11) 



J red 
(12) 



f red 
(13) 



*19 
20 



125 
t26 
27 
28 



36 

f37 
38 



f43 
44 
45 
t46 
47 



53 
54 



57 
58 



60 
t61 



73 
74 
75 



1226 


7813 ± 830 


5 


.40 


± 





57 


2 


.20 


± 


0. 


23 


14.7 ± 3 


,4 


1 


.23 


± 





01 





.87 


± 





18 


1 


.57 


± 





19 





64 


± 





14 


10 


,4 


± 2 


.5 


4 


.3 


± 1 


.3 





.29 


1316 


14660 ± 891 


12 


.59 


± 





.77 


4 


19 


± 





25 


48.5 ± 13 





2 


33 


± 





02 


2 


83 


± 





.74 


3 


04 


± 





21 


1 


14 


± 





15 


35 


.2 


± 9 


.5 


13 


.2 


± 3 


.9 





.27 


1978 


4745 ± 1099 


5 


.16 


± 


1 


.20 


2 


.03 


± 





47 


14 


± 


3 


.8 


1 


15 


± 


o 


02 


o 


84 


± 


o 


12 


1 


28 


± 


o 


.38 


0. 


75 


± 


0. 


27 


8 


.8 


± 2 


.9 


5 


.2 


± 2 


.0 


o 


.37 


881 


2660 ± 129 


2 


.57 


± 





.12 


1 


.03 


± 


0. 


05 


9 


.2 ± 


1 


.9 





.56 


± 





.01 





.49 


± 





.10 





.87 


± 





04 





16 


± 





03 


7 


,7 


± 1 


.6 


1 


.5 


± 


4 





16 


798 


1110 ± 181 


1 


.29 


± 





.21 





.56 


± 





09 


6 


± 


1 


,7 





20 


± 





01 





21 


± 





05 





33 


± 





07 





23 


± 





05 


3 


.5 


± 1 


.1 


2 


,4 


± 


.8 





.41 


763 


4301 ± 1201 


5 


.20 


± 


1 


45 


2 


.26 


± 





63 


18 


.2 ± 


6 


.9 


1 


.11 


± 





.03 





89 


± 





.23 


2 


01 


± 





.51 





24 


± 





37 


16 


.2 


± 5 


.9 


2 


.0 


± 3 








.11 


731 


3246 ± 598 


3 


.86 


± 





71 


1 


.45 


± 


0. 


27 


14 


.3 ± 


4 


.2 





.85 


± 





.02 





86 


± 





.20 





80 


± 





.22 





65 


± 





16 


7 


.9 


± 2 


.8 


6 


.5 


± 2 


.2 





.45 


1535 


388 ± 117 


1 


.09 


± 


o 


33 


o 


.42 


± 


0. 


13 


5 


± 


1 


.8 


o 


.20 


± 


o 


.00 


o 


.25 


± 


o 


.05 


o 


.26 


± 


o 


.10 





16 


± 





07 


3 


.1 


± 1 


.4 


1 


.9 


± 1 





o 


.38 


1903 


803 ± 355 


2 


.70 


± 


1 


.19 


1 


.07 


± 





47 


9 


.6 ± 


4 


.8 





68 


± 





01 





61 


± 





14 





.55 


± 





38 


0. 


52 


± 


0. 


27 


4 


.9 


± 3.6 


4 


,7 


± 2 


.7 





.49 


1632 


984 ± 198 


2 


.82 


± 





.57 


1 


.15 


± 


0. 


23 


10 


.3 ± 


2 


.8 





.65 


± 





.01 





.61 


± 





.11 





.69 


± 





.19 





46 


± 





14 


6 


.2 


± 2 


.0 


4 


,1 


± 1 


4 





.40 


1231 


376 ± 97 


1 


.98 


± 





.51 





.88 


± 





23 


7 


.0 ± 


2 


.4 





32 


± 





01 





28 


± 





06 





50 


± 





18 





38 


± 





13 


4 


.0 


± 1 


.7 


3 


.0 


± 1 


.3 





.43 


2095 


211 ± 34 


1 


.03 


± 





.17 





.88 


± 





14 


4 


± 


1 








16 


± 





00 





14 


± 





03 





52 


± 





11 


0. 


35 


± 


0. 


08 


2 


.4 


± 


.7 


1 


.6 


± 


.5 





.40 


1154 


218 ± 28 


1 


.00 


± 





.13 





.42 


± 





05 


2 


.8 ± 





6 





12 


± 





00 





09 


± 





02 





22 


± 





04 





20 


± 


0. 


03 


1 


.5 


± 


.4 


1 


.3 


± 


.3 





.48 


1062 


178 ± 30 


1 


.06 







18 





.52 


± 





09 


3 


.4 ± 





.8 





.21 


± 





.00 





.16 


± 





03 





.52 


± 





07 


0. 


00 


± 





05 


3 


.4 


± 


.7 





.0 


± 


.3 





.00 


2092 


103 ± 17 





.68 







11 





23 







04 


2 


.2 ± 





.5 





08 


± 





00 





07 


± 





01 





13 


± 





03 


0. 


09 


± 


0. 


02 


1 


.3 


± 


.4 





.9 


± 


.3 





.39 


369 


179 ± 17 


1 


.51 


= 





14 


1 


.23 


= 


0. 


12 


5 


.5 ± 


1 


,4 





28 


± 


o 


01 


o 


25 


± 


o 


06 





.58 


± 


o 


09 





65 


± 





08 


2 


.6 


± 


.7 


2 


.9 


± 


.8 


o 


.53 


759 


200 ± 41 


1 


.34 


= 





27 





.54 


= 


0. 


11 


4 


.4 ± 


1 


3 





.24 


± 





.00 





.20 


± 





05 





35 


± 





.09 





19 


± 


0. 


07 


2 


.9 


± 1 


.0 


1 


6 


± 


.6 





.35 


1692 


139 ± 23 


1 


.09 


= 





18 





.53 


= 


0. 


09 


4 


1 ± 


1 


.0 





16 


± 





01 





16 


± 





03 





34 


± 





07 


0. 


19 


± 


0. 


05 


2 


.6 


± 


.7 


1 


.5 


± 


.5 





.36 


1030 


345 ± 80 


2 


.27 


= 





53 





.97 




0. 


22 


25.9 ± 18 


6 





40 


± 





01 


1 


06 


= 





72 





73 


± 





18 





24 


± 





13 


19.5 ± 14 


.1 


6 


,4 


± 5 


.6 





.25 


2000 


205 ± 28 


2 


.97 


= 





41 


1 


. 12 




0. 


15 


8 


.6 ± 


1 


9 





53 


± 


0- 


01 





.37 


± 





.06 





97 


± 





.13 





15 


± 





08 


7 


.5 


± 1 


.6 


1 


.2 


± 


.7 





.13 


685 


196 ± 60 


1 


.78 


= 





. 54 


o 


.75 


= 


0. 


23 


6 


.3 ± 


2 


.3 





.40 


± 


0- 


01 





36 


= 


0- 


08 





54 


± 


o 


.19 





22 


± 





13 


4 


.5 


± 1 


.8 


1 


.8 


± l 


.2 


o 


.29 


1664 


213 ± 31 


2 


35 


= 





34 





.87 


= 


0. 


13 


8 


.2 ± 


2 


.3 





12 


± 





01 





37 







09 





54 


± 





10 





33 


± 





07 




.1 


± 1 


.5 


3 


,1 


± 1 


.0 





.38 


654 


45 ± 16 





.50 


= 





.18 





.26 


= 


0. 


09 


2 


.0 ± 





.8 





05 


± 





00 





05 


= 





01 





23 


= 





08 


0. 


04 


± 





05 


1 


.7 


± 


.7 





.3 


± 


4 





.15 


944 


72 ± 10 





.87 







12 





.39 




0. 


05 


2 


.5 ± 





.5 





11 


± 


0- 


01 





08 




0- 


01 





25 







.04 





14 


± 





03 


1 


.6 


± 


.4 





.9 


± 


.3 





.36 


1938 


59.8 ± 9.2 





64 


= 





10 





.27 


= 


0. 


04 


2 


.4 ± 





.6 





12 


± 





01 





11 


± 





02 





22 


± 





04 





05 


± 





02 


1 


.9 


± 


.5 





.5 


± 


.2 





.19 


1279 


58 ± 11 


o 


.72 


= 


o 


15 


o 


30 


= 


0. 


06 


2 


.7 ± 





.7 


o 


13 


± 


o 


01 


o 


11 


= 





02 


o 


21 


± 


o 


05 


0. 


09 


± 


0. 


04 


1 


.8 


± 


.6 


o 


.8 


± 


.3 


o 


.30 


1720 


62 ± 13 





.87 







18 





.33 




0. 


07 


3 


.1 ± 


1 


.0 





12 


± 


0- 


00 





10 


± 





.03 





.22 


± 





06 


0. 


12 


± 





04 


2 


.0 


± 


.7 


1 


,1 


± 







.36 


355 


100 ± 31 


1 


.70 


= 





53 


o 


.69 


= 


0. 


21 


6 


.3 ± 


2 


.4 





11 


± 





02 


o 


.11 


± 


o 


03 


o 


53 


± 


o 


.17 





16 


± 





12 


4 


.9 


± 2 


.0 


1 


,4 


± 1 


2 


o 


.23 


1619 


84 ± 19 


1 


.52 







34 





70 




0. 


16 


4 


.9 ± 


1 


.4 





14 


± 





01 





10 


± 





02 





52 


± 





13 


0. 


19 


± 


0. 


09 


3 


.6 


± 1 


.1 


1 


.3 


± 


.7 





.27 


1883 


83 ± 25 


1 


.11 


= 





34 





.42 


= 


0. 


13 




.0 ± 


1 


.9 





.16 


± 





.00 





.16 


± 





.04 





.29 


± 





.10 





13 


± 





07 


3 


.4 


± 1 


.5 


1 


.6 


± 


.9 





.31 


1242 


116 ± 24 


2 


.05 


= 





42 





.88 


= 


0. 


18 


7 


.6 ± 


2 


.2 





31 


± 





01 





26 


± 





05 





.57 


± 





15 





31 


± 


0. 


11 


4 


.9 


± 1 


.6 


2 


.7 


± 1 


1 





.35 


784 


50 ± 14 





.95 







27 





.39 




0. 


11 


3 


.0 ± 


1 


.0 





07 


± 





01 





05 


± 





01 





23 


± 





09 


0. 


17 


± 


0. 


07 


1 


.7 


± 


.7 


1 


.3 


± 


.6 





.44 


1537 


31-4 ± 7-2 





.80 


= 





18 





.33 


— 


0. 


08 


3 


.1 ± 


1 


.0 





00 


± 





00 





00 


± 





00 





29 


± 





07 





04 


± 





03 


2 


.7 


± 


.9 





,4 


± 


.3 





.12 


778 


74 ± 32 


1 


.11 


= 





48 





.72 




0. 


31 


4 


.1 ± 


2 


.0 





14 


± 





.01 





.12 


± 





03 





51 


± 





.25 





21 


± 





18 


2 


.9 


± 1 


.6 


1 


.2 


± 1 


1 





.29 


1321 


31-0 ± 9-0 





.94 


= 





27 





.33 




0. 


10 


4 


.4 ± 


1 


.7 





09 


± 





01 





08 


± 





02 





24 


± 





08 


0. 


09 


± 


0. 


05 


3 


.3 


± 1 


.3 


1 


.2 


± 


.8 





.27 


828 


69-5 ± 9-8 


1 


.52 


= 





21 





.65 


= 


0. 


09 




.1 ± 


1 


.3 





.00 


± 





.00 





.00 


± 





.00 





.50 


± 





.08 





15 


± 





05 


3 


.9 


± 1 


.1 


1 


.2 


± 


4 





.23 


1250 


20-1 ± 7-3 





.52 


= 





19 





.18 


= 


0. 


07 




.5 ± 


3 


.6 





.06 


± 





.01 





.14 


± 





.08 





.17 


± 





06 





01 


± 





03 




.0 


± 3 


.3 





.4 


± 1 


.0 





.06 


1630 


47 ± 11 


1 


.11 







26 





.44 


= 


0. 


10 


4 


.0 ± 


1 


.3 





.12 


± 





.02 





.11 


± 





03 





.26 


± 





.08 


0. 


18 


± 





06 


2 


.4 


± 


.9 


1 


,7 


± 


.7 





41 


1146 


72 ± 12 


1 


93 







32 


1 


.00 


= 


0. 


17 


8 


.3 ± 


2 


.5 





.39 


± 





02 





40 


± 





10 





30 


± 





12 





70 


± 





12 


2 


.5 


± 1 


.1 




.8 


± 1 


.8 





.70 


1025 


141 ± 34 


2 


.08 


= 





50 





.84 


= 





20 


6 


.6 ± 


2 


.1 





.25 


± 





.07 





.19 


± 





08 





.75 


± 





.17 





09 


± 





11 




.9 


± 1 


.8 





.7 


± 


9 





.11 


1303 


72 ± 18 


2 


.10 


= 





53 





.87 


= 


0. 


22 


6 


.9 ± 


2 


.3 





24 


± 





03 





18 


± 





05 





82 


± 





18 


0. 


05 


± 


0. 


11 


6 


.5 


± 2 


.1 





.4 


± 


.9 





.06 


1913 


71 ± 14 


2 


.02 


= 





40 





.90 


= 


0. 


18 


6 


.6 ± 


2 


.1 





31 


± 





01 





23 


± 





06 





74 


± 





15 





16 


± 





09 


5 


.4 


± 1 


.8 


1 


.2 


± 


.8 





.18 


1327 


11 ± 12 





.26 


= 





29 





.14 


± 





15 





.7 ± 





.7 





.07 


± 





00 





.05 


± 





.01 





09 


± 





.12 


0. 


05 


± 





09 





.4 


± 


.6 





.2 


± 


.4 





.36 


1125 


52-3 ± 8-5 


1 


.23 


± 





20 





.88 


± 


0. 


14 


6 


.5 ± 


2 


9 





.16 


± 





01 





.20 


± 





.08 





81 


± 





.13 





07 


± 





06 


6 


.0 


± 2 


.7 





.5 


± 


.5 





.08 


1475 


81 ± 10 


3 


.05 







40 


1 


.38 


± 





18 


10 


.5 ± 


2 


.5 





.33 


± 





.03 





.23 


± 





05 


1 


33 


± 





.17 





06 


± 





06 


10 


.1 


± 2 


.4 





,4 


± 







.04 


1178 


25-3 ± 9-2 


1 


.16 







42 





.46 


± 





17 


3 


.6 ± 


1 


.5 





12 


± 





01 





09 


± 





02 





31 


± 





14 


0. 


15 


± 


0. 


10 


2 


.4 


± 1 


.2 


1 


.2 


± 


.8 





.33 


1283 


58.6 ± 9-3 


2 


.02 







32 





.77 




0. 


12 


6 


.5 ± 


1 


.6 





.22 


± 





03 





.16 


± 





.04 





.59 


± 





.10 





18 


± 





06 


5 


.0 


± 1 


.3 


1 


.5 


± 


6 





.23 


1261 


35-1 ±7-6 


1 


.51 


= 





33 





72 


= 


0. 


16 


7 


2 ± 


2 


9 





29 


± 





02 





28 


± 





10 





60 


± 





14 





12 


± 





08 


6 


.0 


± 2 


.5 


1 


.2 


± 


.9 





.17 


698 


114 ± 12 


3 


.53 


= 





38 


1 


50 


= 


0. 


16 


12 


.0 ± 


3 


.5 





56 


± 





01 





41 


= 





11 


1 


23 


± 





14 


0. 


27 


± 


0. 


07 


9 


.8 


± 2 


.9 


2 


.2 


± 


.8 





.18 


1422 


24.9 ± 6-0 


1 


.62 







39 





80 


= 





19 


6 


.5 ± 


2 


.8 





21 


± 





02 





17 


± 





06 





54 


± 





16 





26 


± 





11 


4 


.4 


± 2 


.0 


2 


,1 


± 1 


.2 





.32 


2048 


17-2 ± 5.4 


1 


.25 







39 





.64 




0. 


20 


5 


9 ± 


2 


8 





19 


± 





02 





.19 


± 





.07 





.49 


± 





.17 





15 


± 





11 


4 


.5 


± 2 


.3 


1 


.3 


± 1 


.1 





.23 


1871 


10-4 ± 5-0 


1 


.24 


= 





60 





49 


= 


0. 


24 


4 


.6 ± 


2 


.5 





21 


± 





07 





17 


= 





07 





48 


± 





21 


0. 


01 


± 


0. 


11 


4 


.5 


± 2 


.3 





.1 


± 1 


.0 





.02 


9 


25-7 ± 6-4 


1 


.57 


= 





39 





.78 


= 


0. 


19 


8 


.3 ± 


2 


.6 





21 


± 





02 





20 


= 





05 





68 


± 





.17 





09 


± 





09 


7 


.3 


± 2 


.4 


1 


.0 


± 1 








12 


575 


18-0 ± 6-1 





.77 







26 





33 




0. 


11 


3 


.6 ± 


1 


,7 





00 


± 





06 


-0 


00 


± 





.06 





.26 


± 





.10 





07 


± 


0. 


06 


2 


.8 


± 1 


.4 





.8 


± 


.7 





21 


1910 


48.7 ± 8.4 


5 


.38 


= 





93 


2 


.27 


= 


0. 


39 


23 


.2 ± 


9 


.8 


1 


09 


± 





04 


1 


12 


± 





43 


1 


47 


± 





32 


0. 


79 


± 


0. 


23 


15 


.1 


± 6 


.7 


8 


.1 


± 3 


.9 





.35 


1049 


8.4 ± 4.5 


1 


.77 


= 





95 





97 


= 





52 


16.2 ± 15 


.1 





05 


± 





03 





09 


= 





08 





.75 


± 





44 





22 


± 





28 


12.4 ± 12.0 


3 


.7 


± 5 


.5 





.23 


856 


43.4 ± 7.9 


4 


.07 


= 





74 


2 


.17 


= 


0. 


39 


19.5 ± 


8 


.3 





55 


± 





03 





.56 


± 





22 


1 


80 


± 





.35 





37 


± 





19 


16 


.2 


± 6 


.9 


3 


.3 


± 2 


1 





.17 


140 


21.3 ± 6.1 


2 


.10 







60 


1 


.05 




0. 


30 


9.2 ± 3 


.8 





.18 


± 





03 





.16 


± 





.06 





.79 


± 





.25 





26 


± 





16 


6 


.9 


± 3 


.1 


2 


.3 


± 1 


6 





.25 


1355 


10.8 ± 5.6 


1 


.07 


= 





.56 





.59 


± 





30 


5 


9 ± 3 


3 





05 


± 





01 





07 


± 





02 





53 


± 





26 


0. 


06 


± 


0. 


15 


5 


.3 


± 2 


.8 





.6 


± 1 


6 





.10 


1087 


66.0 ± 9.5 


5 


.07 


= 





73 


2 


31 


= 


0. 


33 


20.0 ± 7 


.1 





.87 


± 





.03 





.73 


± 





.24 


2 


17 


± 





.31 





14 


± 





12 


18 


.8 


± 6 


.7 


1 


.2 


± 1 


1 





06 


1297 


4 ± 11 





.42 


= 


1 


00 





.16 


- 


0. 


37 


1.2 ± 3 


.0 





.15 


± 





02 





.12 


± 





.02 





.11 


± 





.30 





05 


± 





22 





.9 


± 2 


,4 





.4 


± 1 


.7 





.31 


1861 


37.6 ± 7.4 


3 


.44 


= 





.68 


1 


.73 


± 





34 


13.0 ± 5 


.2 





69 


± 





07 





49 


± 





18 


1 


26 


± 





29 


0. 


47 


± 





18 


9 


.5 


± 3 


.9 


3 


.5 


± 1 


.8 





.27 


543 


18.1 ± 5.5 


1 


.97 


± 





.60 





.96 


± 





29 


8 


.3 ± 3 


.5 





15 


± 





03 





14 


± 





05 





96 


± 





.27 


0. 


00 


± 


0. 


11 


8 


.3 


± 3.4 





.0 


± 


9 





.00 


1431 


60.6 ± 9.3 


6 


.72 


± 


1 


.03 


3 


.03 


± 





46 


27.5 ± 10 


.2 





.78 


± 





04 





.70 


± 





.24 


2 


.57 


± 





.41 





45 


± 





21 


23 


.4 


± 8 


.7 


4 


.1 


± 2 


.4 





15 


1528 


40.7 ± 7.6 


5 


.57 


± 


1 


04 


2 


.46 


± 





46 


24 


9 ± 


8 


.7 





66 


± 





04 





69 


± 





21 


2 


24 


± 





42 





22 


± 





19 


22 


.7 


± 7 


.9 


2 


.2 


± 2 








.09 


1695 


14.4 ± 5.7 


1 


.45 


± 





.58 





.68 


± 





27 


8 


. 5 ± 5 


.2 





09 


± 





03 





10 


± 





06 





48 


± 





22 





20 


± 


0. 


15 


6 


.0 


± 3 


.9 


2 


.5 


± 2 


.2 





.29 


1833 


18.1 ± 5.5 


2 


.55 


± 





77 


1 


.26 


± 


0. 


38 


21.1 ± 13 


.5 





.15 


± 





.04 





.25 


± 





.16 





.87 


± 





.32 





39 


± 





22 


14 


.6 


± 9 


.8 


6 


.5 


± 5 


.2 





.31 


437 


42.1 ± 7.9 


3 


.13 


± 





.59 


1 


.40 


± 





26 


15.0 ± 7 


.2 





46 


± 





04 





34 


± 





15 


1 


32 


± 





24 


0. 


08 


± 


0. 


10 


14 


.1 


± 6 


.8 





.9 


± 1 


.2 





.06 


2019 


23.9 ±6.1 


2 


.72 


± 





.69 


1 


.28 


± 





33 


23.5 ± 17 


.8 





.72 


± 





05 


1 


19 


± 





.85 


1 


.28 


± 





31 


0. 


00 


± 


0. 


11 


23.5 ± 17.6 





.0 


± 2 


.1 





.00 


33 


2.2 ± 4.2 





.61 


± 


1 


17 





.34 


± 





66 


5.1 ± 10 


.1 


-0 


13 


± 





07 


-0 


21 


± 





16 





22 


± 





53 


0. 


12 


± 





38 


3 


.3 


± 8 


.1 


1 


.8 


± 5 


7 





.35 


200 


15.5 ± 5.8 


2 


.21 


± 





.83 


1 


.18 


± 





44 


11.5 ± 6 


.0 





.05 


± 





.05 





.05 


± 





.06 


1 


18 


± 





.40 


0. 


00 


± 





19 


11.5 ± 5.8 





.0 


± 1 


.8 





.00 


571 


10.9 ± 5.6 


1 


.29 


± 





.66 





.53 


± 





27 


11.2 ± 9 


.4 





14 


± 





07 





17 


± 





20 





53 


± 





24 


0. 


00 


± 


0. 


13 


11.2 ± 9.0 





.0 


± 2 


7 





.00 


21 


20.7 ± 6.5 


3 


.85 


± 


1 


21 


2 


.32 


± 


0. 


73 


37.0 ± 27 


.2 





.67 


± 





.08 





97 


± 





.65 


2 


32 


± 





.66 





00 


± 





30 


37.0 ± 26.7 





.0 


± 4 


.8 





.00 


1488 


7.3 ± 4.2 


1 


.42 


± 





.82 





.85 


± 


0. 


49 


17. £ 


1 ± 19 


.5 





.03 


± 





04 





.08 


± 





12 





84 


± 





.44 


0. 


01 


± 





22 


17.7 ± 18 


.8 





.2 


± 4 


6 





.01 


1779 


2.1 ± 3.7 





.37 


± 





.64 





.23 


± 





41 


2 


.1 ± 3 


.7 


-0 


.11 


± 





.05 


-0 


.13 


± 





08 





.23 


± 





35 





00 


± 





22 


2 


.1 ± 3.2 





.0 


± 2 


.0 





00 


1895 


6.3 ± 4.3 


1 


.45 


± 





99 





.76 


± 





52 


8 


.3 ± 7 


.2 





07 


± 





06 





05 


± 





07 





65 


± 





44 


0. 


11 


± 





27 


7.1 ± 6 


.1 


1 


.3 


± 3 


.1 





.14 


1499 


19.0 ± 6.6 


4 


.65 


± 


1 


62 


3 


.14 


± 


1. 


09 


133.8 ± 185 


1 





98 


± 





14 


4 


.07 


± 


5 


47 


2 


67 


± 





94 





47 


± 





55 


113.7 ± 157.5 


20.] 


. ± 35 


.7 





.15 



o 

p 

cro 
CD 



TABLE 2 — Continued 



o. 

(i) 


( 2 ) 


N GC 
(3) 




S N 
(4) 








S N,z 
(5) 




(6) 






s L 

(7) 








(S) 






S N,z,blu 
(9) 






N,z,re 
(10) 


d 


1 blue 
(11) 




^Ved 
(12) 




J red 
(13) 


78 


1545 


54.2 zt 8.8 


9 


37 


± 


1 


52 


4 


.34 


zb 


.71 


38.5 zb 13 


o 


1 


.02 


zb 


Q 


07 


o 


.76 zb 


o 


23 


3 


91 zb 





64 


o 


.43 


zb 


.29 


34.6 zb 11 


8 


3.8 


± 2 


8 


0. 10 


f79 








12 


± 








.34 


zb 










39 


zb 








31 zb 






— 


21 zb 









.13 


zb 








— 1.3 


± 1 






80 


1857 


10.8 ± 5.9 


2 


.39 


± 


1 


.31 


1 


.28 


zb 


.70 


14.8 zt 8 


Q 


o 


32 


zb 


o 


09 


o 


29 i 


o 


12 





78 ± 


o 


56 


o 


.50 


zb 


.42 


9.1 zb 6 


7 


5.8 


± 5 





0.39 


81 


1075 


16.5 ± 5.2 


3 


19 


± 


1 


01 


1 


.58 


zb 


.50 


14.8 zb 7 


g 


o 


.33 


zb 


o 


.06 


o 


31 zb 


o 


14 


1 


58 ± 


o 


. 46 


o 


.00 


zb 


. 18 


14.8 zb 7 


Q 


0.0 


± 1 


7 


0.00 


82 


1948 


2.5 zt 3.0 


o 


.9-1 


± 


1 


. 13 


o 


.48 


zb 


.58 


5.0 zb 6 


1 


— 


28 


zb 


o 


17 


— 


22 zt 


o 


12 


o 


48 ± 


o 


50 


o 


.00 


zb 


.29 


5.0 ± 5- 


3 


0.0 


± 3 


.0 


0.00 


83 


1627 


3.6 zt 3.7 


o 


.94 


± 


o 


97 


o 


.40 


zb 


.41 


3 .5 zb 3 


8 


1 


21 


zb 


o 


13 


1 


05 zb 


o 


34 


o 


40 zb 


o 


36 


o 


.00 


zb 


,21 


3.5 zt 3 


3 


0.0 


± 1 


8 


0.00 


84 


1440 


26.7 zt 6.8 


4 


81 


= 


1 


22 


2 


.30 


zb 


.58 


22.5 zt 10. 


1 


o 


.29 


zb 


o 


05 


o 


-31 ± 


o 


. 12 


1 


93 zb 


o 


.51 


o 


.37 


zb 


.28 


18.9 zt 8 


Q 


3.6 


± 3 


1 


0. 16 






28.7 zt 6.7 




.39 


± 








.73 


± 1 










.45 


zb 














4 


64 ± 


1 






.09 


zb 




40. 5 zt 19 


* 


0.8 


± 3 






86 


2050 






63 


± 




23 




.27 


± 


60 


31 5 zb 32 






.21 


zb 




06 


j? 


43 zb 




41 


1 


.27 zb 


o 


54 




.00 


zb 


25 


^ 1 'o ^±2 




0.0 


± 6 




00 


87 


1993 






.48 


± 









.25 


zb 


33 


2 1 ± 2 







07 


zb 




08 





29 zb 





13 


-0 


03 ± 





30 





.22 


zb 


99 






-1.9 


± 2 


6 


00 


88 


751 


9.2 ± 4.3 


i 


49 


± 


o 


70 


o 


.62 


zb 


.29 


6.5 zt 3 


g 


o 


.13 


± 


Q 


04 


Q 


. 11 zb 


o 


06 





62 ± 





27 


o 


.00 


zb 11 


12 


6.5 ± 3 


Q 


0.0 


± 1 


3 


0. 00 


89 


1828 


20.4 ± 5.8 


3 


94 


± 


1 


.12 


2 


.10 


± 


.60 


16.3 zb 7 


7 


o 


95 


± 


o 


08 


o 


78 ± 


o 


30 


2 


10 =b 





56 


o 


.00 


± 


.22 


16.3 zt 7. 


5 


0.0 


± 1 


.7 


0.00 


90 


538 


0.6 ± 3.5 





19 


= 


1 


.11 


o 


.09 


± 


.51 


0.9 zb 5 


5 


— 


31 


± 


o 


22 


— 


25 ± 


o 


19 





09 zb 





42 


o 


.00 


zb 


.29 


0.9 zt 4 


6 


0.0 


± 3 


1 


0.00 


91 


1407 


49.7 ± 8.6 


10 


18 


= 
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Running number, sorted by increasing magnitude 

2 

Number in Virgo Cluster Catalog 

3 

Total number of globular clusters 

4 

Specific frequency 

Specific frequency in z bandpass 
6 Nq C normalized to stellar mass of 10 9 Mq 

7 

Percentage of galaxy z luminosity in GCs 
Percentage of galaxy stellar mass in GCs 

9 

Specific frequency in z for blue GCs 
Specific frequency in z for red GCs 

11 T for blue GCs 

12 T for red GCs 

13 

Fraction of red GCs 
^ Quantities calculated using an R.= 70 // aperture around galaxy 

^ Quantities possibly affected by nearby dE VCC 1941, which has been masked 
* Mass-to- light ratio (i.e. optical-IR colors) suspect. 



